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Preface 


The  purpose  of  this  thesis  was  to  derive  a  set  of  equations 
that  exactly  determine  the  scattered  EM  fields  from  a  cavity  in 
an  infinite  conducting  plane.  The  chosen  approach  was  to  use  a 
set  of  integral  equations  for  the  field  components  inside  the 
cavity  and  a  set  of  integral  equations  for  the  field  components 
outside  the  cavity  and  couple  them  together  at  the  junction.  A 
Green’s  function  with  Neumann  boundary  conditions  was  used 
outside  the  cavity.  The  thesis  contains  derivations  of  all  the 
major  equations  used. 
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Abstract 


The  purpose  of  this  thesis  is  to  determine  the  scattered 
electromagnetic  fields  and  radar  cross  section  from  a  two- 
dimensional  cavity  in  a  perfectly  electric  conducting  infinite 
plane.  This  is  accomplished  by  deriving  a  coupled  set  of 
Fredholm  integral  equations  of  the  second  kind.  A  set  of 
integral  equations  outside  the  cavity  and  a  set  of  integral 
equations  inside  the  cavity  are  coupled  together  at  the 
interface.  The  Fredholm  integral  equations  of  the  second  kind 
for  the  outside  of  the  cavity  use  a  Green’s  function  with  Neumann 
boundary  conditions  to  avoid  an  integration  over  the  infinite 
plane  for  a  transverse  electric  incident  plane  wave.  An  example 
problem  is  introduced  and  numerically  solved  to  test  the 
application  of  the  newly  derived  equations. 


Coup  1 ed  Integral  Equation  So  1 ut ions  f or  Two-Diaens ional 
Bistatic  TE  Scatter  From  a  Conducting 
Cavity-Backed  Inf  inite  Plane 

I_.  Introduction 

Electromagnetic  (EH)  waves  exhibit  the  phenomena  of 
reflection,  transmission,  absorption,  interference, 
diffraction,  and  refraction  (13,11,4).  In  free  space, 
absorption,  diffraction,  and  refraction  may  be  ignored  (14). 
With  this  single  and  precise  assumption,  EM  energy  (such  as  a 
radar  pulse)  propagates  through  space  until  it  impinges  upon  an 
obstruction  (such  as  an  aircraft)  in  its  path.  The  obstruction 
scatters  the  electromagnetic  wave  in  many  directions  and  the 
energy  associated  with  the  scattered  portion  of  the  EM  wave 
then  propagates  in  new  directions.  Though  the  energy  is 
normally  scattered  in  all  directions,  it  is  desirable  to  know 
the  amount  of  energy  scattered  in  various  directions. 

Solutions  to  this  problem  in  the  form  of  the  scattered  electric 
and  magnetic  fields  are  of  interest  to  many  sectors  of  today’s 
society,  particularly  to  the  military  (6). 

It  is  not  a  trivial  problem  to  determine  the  exact 
solution  for  the  scattered  EH  fields.  The  size  and  shape  of 
the  obstruction  must  be  taken  into  account  when  solving  for  the 
scattered  fields,  but  these  are  only  two  of  the  possible 
variables.  The  object’s  conduct ivity,  the  incident  wave’s 
frequency,  and  the  permeability  and  permittivity  of  the  space 


surrounding  the  object  all  contribute  to  the  solution  of  the 
scattered  fields  (6). 

Closed  form  solutions  of  the  EM  scattering  problem  can  be 
generated  for  only  a  very  limited  class  of  problems  (12,  14). 

For  general  scattering  objects,  however,  an  approximate 
solution  is  needed  to  predict  the  scattered  fields.  With  high 
speed  computer  systems  (or  low  cost  personal  type  computers ) 
good  numerical  approximations  are  obtainable  by  solving  the 
full  scattered  field  equations  (12). 

The  Prob 1 em 

Modern  jet  aircraft  use  air  intake  and  exhaust  structures 
that  can  appear  (to  the  EM  wave)  to  be  open  ended  wave  guides 
or  apertures.  If  one  wishes  to  reduce  the  energy  bacKscattered 
by  these  apertures,  it  is  necessary  to  determine  both  the 
amount  of  EM  energy  scattered  by  them  and  the  directions  into 
which  the  energy  is  scattered.  This  thesis  makes  another  step 
toward  unders tand  ing  this  three-dimensional  problem  by 
determining  the  solution  to  a  corresponding,  but  simpler,  two- 
dimensional  problem  (26).  One  benefit  from  solving  this 
simpler  problem  is  the  two-dimensional  problem  can  be  used  to 
approximate  the  scattering  from  gaps  and  cracks  in  panels  that 
comprise  reflector  antennas  and  radar  targets. 

The  associated  two-dimensional  model  of  the  problem  that 
will  be  solved  in  this  thesis  is  the  determination  of  the 
bistatic  EM  fields  scattered  from  an  infinite  plane  containing 
a  cavity  (26).  For  example,  if  in  three  dimensions  the 


2 


infinite  plane  is  the  xz  plane  and  the  coordinate  axes  are 
placed  as  in  Figure  i,  then  the  cross-sectional  shape  of  the 
cavity  is  constant  along  the  z-axis.  This  is  the  specific 
problem  ultimately  considered  in  this  thesis. 

Integral  equations  can  be  derived  to  describe  the 
electromagnetic  fields  scattered  in  a  problem  of  this  type. 
These  equations  are  developed  in  this  thesis.  Once  derived, 
computer  software  will  be  described  which  numerically 
implements  this  formulation.  In  particular,  predictions  for 
two-dimensional  scattering  from  an  infinite  plane  containing  a 
rectangular  shaped  cavity  will  be  obtained. 

Scope  of  the  Prob 1 em 

It  is  desired  to  determine  the  EM  fields  scattered  from  a 
perfectly  electric  conducting  infinite  plane  containing  a 
cavity  of  arbitrary  cross  section  (i.e.  a  generic  cavity). 
Surface  integral  equations  have  been  derived  by  Stratton, 
Poggio  and  Miller  and  others  (21,  15,  23,  24)  that  provide  for 

the  solution  of  scattered  fields  from  objects  of  finite  size 
using  a  free  space  Green's  function  as  the  kernel.  To  avoid 
integrations  over  the  unbounded  infinite  plane,  Green’s 
functions  which  satisfy  Dirichlet  boundary  conditions  (the 
Gre.'n’s  function  vanishes  at  the  boundary)  or  Neumann  boundary 
conditions  (the  normal  derivative  of  the  Green’s  function 
vanishes  at  the  boundary)  can  be  employed. 

The  integral  equations  using  a  Green’s  function  with 
Neumann  boundary  equations  will  be  derived  and  used  in  this 
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As sump t ions 


The  problem  is  simplified  by  assuming  that  the  infinite 
plane  is  a  perfect  electric  conductor  (i.e.  the  conductivity  of 
the  plane  is  infinite).  The  plane  is  suspended  in  free  space 
where  the  incident  wave  will  not  be  impeded  (i.e.  the  relative 
permeability,  pr,  and  the  relative  permittivity,  er  are 
both  equal  to  one,  pr  =  eP  =  1).  Though  a  perfect 
conductor  of  this  type  does  not  actually  exist,  the 
mathematical  assumption  of  a  PEC  (perfect  electric  conductor) 
is  a  good  approximation  and  the  amount  of  work  involved  in 
assuming  the  material  is  not  perfect  is  considerably  more 
difficult  ( 6  ) . 

The  assumption  that  the  relative  permeability  and 
permittivity  are  exactly  equal  to  one  means  that  the  medium 
through  which  the  EM  wave  propagates  does  not  impede  the  wave 
(14).  In  actuality,  €P  for  air  is  approximately  1.0006  at 
atmospheric  pressure  and  24°C  (11:58)  and  pr  for  air  is 
approximately  1.0000004  at  the  same  conditions  (11:216).  Since 
both  €r  and  pr  are  functions  of  temperature,  pressure,  and 
material,  they  will  be  functions  of  position  in  space.  It  is 
not  possible  to  determine  the  exact  permittivity  and 
permeability  for  all  places  at  all  times  but  it  is  known  that 
for  air,  they  will  be  very  close  to  one  (14). 

One  final  assumption  for  solving  this  problem  is  that  the 
source  of  EM  waves  is  far  enough  away  that  the  wave  front  is 
planar  as  it  impinges  upon  the  plane. 
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Approach 


The  general  geometry  initially  considered  in  this  problem 
is  a  three  dimensional  generic  cavity  in  a  PEC  (perfectly 
electrically  conducting)  infinite  plane  (see  Figure  1).  Within 
the  thesis  the  specific  geometry  of  a  two-dimensional 
rectangular  channel  is  introduced  as  the  cavity,  the  integral 
equations  inside  and  outside  of  the  cavity  are  then  solved  to 
produce  the  fields  scattered  from  the  channel.  For  purposes  of 
derivation,  the  cavity  will  be  considered  generic.  The  chosen 
approach  to  solving  this  problem  is  to  separate  the  upper  half 
space  from  the  cavity  and  allow  the  two  regions  to  be  coupled 
at  their  common  boundary,  hereafter  referred  to  as  the 
"aperture".  This  method  was  introduced  by  Harrington  and  Mautz 
in  1976  and  lends  itself  well  to  this  type  of  problem  (5). 

By  considering  the  scatterer  to  be  a  PEC  plane  with  a 
cavity,  two  sets  of  equations  need  to  be  derived:  one  set  for 

the  half  space  above  the  PEC  plane  and  one  set  for  the  cavity. 
The  first  step  used  to  solve  this  problem  was  to  derive  from 
Maxwell’s  Equations  a  set  of  integral  equations  needed  to 
determine  currents  (field  components)  in  the  aperture  and  along 
the  cavity  walls.  From  the  currents,  the  scattered  fields  are 
obtained  by  an  integration.  An  integral  equation  formulation 
was  used  because  it  reduces,  by  one  dimension,  the 
dimensionality  of  the  problem  and  represents  an  exact  solution. 
Other  theoretical  approaches  such  as  physical  optics, 
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geometrical  optics,  and  the  geometrical  theory  of  diffraction, 
etc.,  require  major  analytic  approx  imat  i  ons  (6,7,18,19). 

Once  the  integral  equations  are  derived,  a  computer 
program  is  introduced  that  solves  the  equations  for  the  unknown 
magnetic  surface  current  density,  M,  from  various  aspect  angles 
and  scattered  angles.  With  the  magnetic  surface  current 
density  in  hand,  the  scattered  fields  are  then  determined  by  an 
integration  of  M.  An  equation  for  the  RCS  (radar  cross 
section)  is  then  derived  and  evaluated  to  determine  the  RCS  of 
the  channe 1 . 

With  the  integral  equation  formulation  and  the  computer 
program,  data  is  generated  and  plotted  to  show  the  RCS  as  a 
function  of  several  variables  for  a  representative  chosen 
geometry:  that  of  a  two-dimensional  rectangular  cavity.  The 

AFIT  computer  facilities  were  used  to  run  the  computer  program. 
A  personal  computer  was  used  to  plot  the  data  generated  from 
the  computer  program. 

Before  preceding  to  the  development  of  the  equrtions,  a 
summarization  of  the  current  knowledge  is  necessary. 

Summary  of  Current  Knowl edge 

Integral  equations  for  finding  the  electric  and  magnetic 
fields  scattered  from  a  three-d imens iona 1  object  have  been 
derived  and  used  for  several  problems.  Two  of  ti.t  most  notable 
derivations  are  by  J.A.  Stratton  (20)  and  Poggio  and  Miller 
(15).  These  equations  are  normally  called  the  Electric  Field 
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Integral  Equation  (EFIE)  and  the  Magnetic  Field  Integral 
Equation  ( MFIE )  (19). 

In  three  dimensions,  both  the  EFIE  and  the  MFIE  contain  a 
singular  kernel  in  convergent  integrals;  however  the  MFIE  is 
normally  more  numerically  stable  and  slightly  less  complicated 
to  implement  (22).  For  open  conductors  such  as  an  antenna 
aperture,  the  MFIE  becomes  degenerate  and  the  EFIE  is  normally 
used.  In  the  early  1980's  the  EFIE  and  the  MFIE  were  combined 
or  augmented  to  eliminate  spurious  resonances  which  can  pollute 
the  solution  of  an  exterior  scattering  problem  (24,  23).  These 
integral  equations  may  be  used  to  determine  the  scattered  EM 
fields  for  many  types  of  problems. 

It  has  been  possible  to  determine  the  scattered  fields 
from  some  of  the  simplest  geometries  since  just  after  the 
introduction  of  Maxwell's  electromagnetic  field  equations  (11). 
Only  recently,  high  speed  computers  have  allowed  the  use  of  the 
MFIE  and  EFIE  to  determine  the  fields  scattered  from  more 
complicated  geometries  (12:2).  Several  computer  codes  have 
been  written  using  the  MFIE  or  EFIE  to  help  determine  the 
fields  scattered  from  some  specific  geometries  (6).  Using 
these  codes,  the  scattered  fields  from  several  geometries  have 
been  determined.  Most  of  these  are  simple  three  dimensional 
bodies.  The  open  literature  does  not  provide  any  solutions  to 
most  of  the  more  complex  shapes.  Some  representations  of 
scattered  fields  that  have  been  obtained  using  surface  integral 
equations  are: 
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1)  The  perfectly  conducting  cube  for  monostatic 
scattering  using  both  the  MFIE  and  the  EFIE  (25). 

\ 

2)  An  infinitely  long,  perfectly  conducting  circular 
cylinder  for  monostatic  scattering  using  both  the  EFIE  and 
the  MFIE  (21  )  . 

3)  An  infinite  perfectly  conducting  half  plane  with  short 
"bumps”  (16)  for  monostatic  scattering  and  using  only  the 
EFIE. 

4)  Various  generic  targets  within  stratified  media  (as 
opposed  to  free  space)  using  only  the  MFIE  (17). 

5)  Perfectly  conducting  bodies  of  revolution  using  both 
the  EFIE  and  the  MFIE  in  the  monostatic  case  (8). 

In  order  to  use  the  MFIE  or  EFIE  in  our  stated  problem  of 

interest,  either  an  integration  on  an  unbounded  domain  is 

necessary  or  else  an  approximation  must  be  made  by  truncating 

the  plane  at  some  finite  distance  from  the  cavity.  The 

necessary  approximation  is  not  considered  valid  in  an  "exact" 

solution.  Therefore,  a  set  of  equations  are  derived  herein 

that  avoids  the  unbounded  domain  without  the  use  of  an 

approx imat ion . 

Preview 

A  coupled  set  of  Fredholm  integral  equations  of  the  second 
kind  are  derived  in  Chapter  2.  Once  obtained,  these  integral 
equations  will  be  simplified  in  Chapter  3  by  reducing  them  to 
two-dimensional  space  and  providing  a  specific  incident 
field.  Chapter  4  provides  a  way  to  solve  the  coupled  set  of 
Fredholm  integral  equations  numerically;  the  Method  of  Moments 
is  used.  In  Chapter  5,  a  specific  cavity  is  introduced  to  test 
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the  equations.  The  specific  cavity  is  a  two-dimensional 
rectangular  channel.  A  convenient  representation  of  the 
scattered  fields  is  the  RCS  (radar  cross  section)  and  an 
equation  for  RCS  is  derived  in  Chapter  6.  Chapter  7  provides  a 
quick  description  of  the  computer  program  used  to  approximate 
the  equations  derived  in  Chapters  2  and  3  to  determine  the  RCS 
derived  in  Chapter  6.  Chapter  8  gives  the  results  of  the 
tested  cavity  and  Chapter  9  gives  the  overall  conclusions  and 
recommendations  for  future  work.  Three  appendices  are  included 
at  the  end  of  the  thesis.  Appendix  A  contains  the  computer 
program  written  to  numerically  solve  the  integral  equations  and 
eventually  determine  the  RCS  from  the  test  (rectangular) 
cavity.  Appendix  B  consists  of  some  of  the  output  from  running 
the  computer  program.  A  description  of  a  model  designed  and 
built  to  further  test  the  equations  by  making  laboratory 
measurements  of  the  scattered  fields  is  in  Appendix  C. 

The  first  major  step,  however  is  to  derive  the  coupled  set 
of  Fredholm  integral  equations  needed  to  solve  for  the  surface 
currents  and  charges,  and  this  is  the  subject  of  the  next 
chapter . 
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1 1 ■  Per ivat ion  of  the  Integral  Equations 


The  Upper  Ha  1 f  Space  Region 

The  propagation  of  electromagnetic  energy  and  the 
fundamental  laws  of  electromagnetics  are  governed  by  Maxwell’s 
equations  (4,11,14,20).  To  determine  the  magnetic  and  electric 
fields  at  any  point  x  =  (x,  y,  2)  at  a  time,  t,  one  may  begin 
with  Maxwell’s  equations  and  derive  a  more  useful  set  of 
equations  for  the  particular  problem  being  considered  (20). 

For  chis  reason,  the  integral  equations  used  to  solve  the 
stated  problem  will  be  derived  from  Maxwell’s  equations. 

In  time  harmonic  form  with  exp(-iwt)  time  dependence, 
Maxwell’s  equations  are  (11:376): 


V  x  E  =  lUJpH  -  I 

(1  ) 

V  x  H  =  -iweE  +  J 

(2  ) 

V  •  E  =  p/e 

(3  ) 

V  •  H  =  pm/p 

(4) 

where 


E  =  Electric  Field  Vector 
H  =  Magnetic  Field  Vector 
I  =  Magnetic  Current  Volume  Density 
J  =  Electric  Current  Volume  Density 
u>  =  Operating  Angular  Frequency 
p  =  permeability  of  the  medium 
e  =  permittivity  of  the  medium 
p  =  Electric  Charge  Volume  Density 
pm  =  Magnetic  Charge  Volume  Density 


All  of  the  fields  (E  and  H),  currents  (I  and  J),  and 
charges  (p  and  pm )  in  Eqs  (1)  through  (4)  are  functions 
of  position  (i.e.  E  =  E(x,  y,  z  )  =E(x  ),  H  =  H(x,  y,  z  )  =H(x  ),  etc.)  For 
simplicity,  the  position  variables  will  not  be  indicated  until 
it  is  necessary  to  distinguish  between  separate  positions  in 
space . 

At  this  point  it  should  be  noted  that  the  magnetic 
current,  I,  and  the  magnetic  charge,  pm,  are  fictitious 
quantities  that  can  exist  mathematically,  but  have  not  yet  been 
physically  measured.  The  electric  and  magnetic  fields  in  Eqs 
(1)  through  (4)  are  considered  generated  from  the  sources  J,  I, 
pm,  and  p.  The  continuity  equations  relating  currents 
and  charges  of  Eqs  (1)  through  (4)  are  (20:464): 

V  •  J  r  iwp  ( 5  ) 

and 

V  •  I  =  iiupm  (6) 

Taking  the  curl  of  Eq  (1)  yields 

V  x  V  x  E  =  iu>p(V  x  H)  -  V  x  I  (7) 

Using  Eq  (2)  in  Eq  (7)  and  simplifying  the  result  gives  an 

equation  for  the  electric  field  without  the  magnetic  field. 

VxVxE-K2E  =  iu>p  J  -  V  x  I  (0) 

where 

k2  =  u)2ey 
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Similarly  for  the  magnetic  field, 


VxVxH-k2H  =  iu>e  l  -  V  x  J  (9) 

An  application  of  a  vector  Green’s  identity  with  the 
proper  Kernel  (or  Green’s  Function)  will  yield  two  integral 
representations  of  the  electric  and  magnetic  fields  in  terms  of 
volume  and  surface  currents  and  charge  densities.  From  these, 
integral  equations  are  derived  to  obtain  the  surface  current 
and  charge  densities.  The  vector  analog  to  the  scaler  Green’s 
second  identity  is  (20:464): 


(  Q  •  VxVxP-P'VxVxQ  )  d  v 


o 

V 


(PxVxQ-QxVxP  )n  ds 


(10  ) 


o  u 
S 


where  n  is  the  unit  normal  pointing  out  of  a  regular  region,  V, 
bounded  by  the  surface,  S.  Eq  (10)  is  valid  for  any  vector 
functions  of  position  P  and  Q  having  the  proper  regularity. 
Both  P  and  Q  need  to  be  twice  continuously  differentiable  on 
the  surface  and  within  the  volume  being  considered  (15:161). 

For  the  vector  P,  either  E  or  H  is  chosen.  Choosing  the 
most  advantageous  vector  Q  depends  upon  the  geometry  being 
considered.  A  generic  closed  surface  in  free  space  lends 
itself  well  to  the  free  space  Green’s  function  for  Q. 
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That  is: 


where 


Q  =  5  (x;  x  *  )  a 


(11  ) 


5(x;  x»  )  = 


r  = 


exp ( iKr ) 
r 

x  -  x  ’ 


(12  ) 


=  £(x-x’  )2  +  (y-y '  )2  +  (z-2  '  )2  j 


1/2 


(13  ) 


a  =  arbitrary  unit  vector 


The  unit  vector,  a,  is  arbitrary  for  two  reasons:  the  free 

space  Green’s  function  is  omnidirectional  and  the  vector  is 
common  to  all  of  the  elements  in  Eq  (10)  (20:251-252).  For  the 

half  space  above  the  PEC  plane,  the  vector  Q  is  chosen  to  be 
the  half  space  Green’s  function: 


3  =  $(x;x’)  a  +  $(x;x’)  a  =  G(x;x’)  a  (14) 

where 

=  (x ’ ,  -y  ’ ,  z  '  )  (15) 

and  0(x;x’)  is  as  defined  in  Eq  (12). 


Using  Eq  (14)  for  Q  and  E(x)  for  P  in  Eq  (10)  yields: 


G  (x;  x  ’  )a  •  V'  x  V'  x  E ( x »  ) 

V  -  E(x')  •  V*  x  V*  x  G  (x;  x  * ) a  j  dv’ 

=  JJ  [e  (x '  )  x  V*  x  G(x;x*)a  -  G(x;x’)a  x  V’  x  E(x’)  j  •  n (x  ’  )  ds 


(16) 


in  which  V1  represents  differentiation  with  respect  to  x’. 

From  Maxwell’s  equations,  V’x  E  is  known  as  a  function  of 
H  and  I  and  Eq  (8)  gives  V’x  V’x  E;  the  only  other  expressions 
that  need  to  be  analyzed  are  V'x  G(x;x’)a  and  V’x  V’  x  G(x;x’)a. 
From  basic  vector  identities: 


V’x  [G  (x;  x  ’  )  a]  =  V’[G(x;x’)]  x  a  (17) 

and 

V’x  V’x  G  (x;  x  ’  )  a  =  V*  [V’  •  (G  (x;  x  ’  )a  )]  -  V’2  [G(x;x’)a]  (18) 

Dealing  strictly  with  Eq  (17)  and  the  integrand  in  the 
left  side  of  Eq  (16)  (recall  that  a  is  a  constant  unit  vector): 

[(E  x  V’  x  (G  a ) ] *  n  =  a-[(n  x  E)  x  V’G)  (19) 

and 


[G(x;x’)a  x  V’x  E(x’)]-n  =  a1  [V’x  E(x’)  x  n)  G(x;x’)  (20) 
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Using  Eq  (i)  in  Eq  (20)  for  V’  x  E(x’)  and  simplifying  the 
result  yields: 

-G(x;x’)a  x  7,x  £(x')-n  =  a-iu>Mtn  x  H(x’)]  G(x;x’) 

-a1  [n  x  I  (x*  )]  G(x;  x*  )  (21  ) 

To  further  evaluate  Eqs  (18)  and  (19),  it  is  necessary  to 
determine  V'G(x;x*)  and  V,2G(x;x*)* 


V'G(x;  x*  )  = 


3G(x;  x*  ) 


3x 


3G (x;  x  '  )  3G  (x;  x  •  ) 

A  "  ■”  A  —  — 

x  +  -  y  ♦ 


3y* 


3z 


(22) 


d  [exp ( ikr ) ] /r 


dr 


3r 


A 

X 


3x 


3r 


ay 


y  + 


3r 


3z 


d  [exp ( ikr ) ] /r 


dr 


3r 


3r 


A 

X  + 


Y  * 


3x 


3y» 


3r 


3z 


(23  ) 


Breaking  Eq  (23)  into  its  prominent  pieces  yields: 


d [exp ( ikr ) ] /r 


dr 


d [exp ( ikr  ) ) /r 

M 

dr 


[“ 

[« 


0(x;  x*  ) 


5  (x;  x‘  ) 


(24) 


(25) 
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and  continuing, 


3r 

(x-x 

’  ) 

(26) 

3x  * 

r 

3r 

(  2-2  ’ 

) 

( 27  ) 

32  ’ 

r 

where,  in  Eqs  (26) 

and 

(27),  r  may  replaced  everywhere 

by  r . 

Continuing, 

dr 

(y-y’ ) 

(28) 

ay’ 

r 

and 

3r 

(y+y* 

) 

(29) 

ay’ 

r 

Thus,  combining 

Eqs  (24)  through  (29)  into  Eq  (23)  yields: 

V’G(x;  x’  )  =  £ik 

-  ■ 

-H 

0 (x;  x 1 ) 

-  r 

r 

*  [ 

ik 

-  — 

1  -I  $(x;x)  r 

.  J  ?  [2  ( y+  y  ’ )  y  -  r  J 

(30) 

where 

r  =  x  - 

x  ’ 

=  (x 

-x')x  +  (y-y* )y  ♦  (2-2’  )2 

m 

r  =  x  - 

x  ’ 

=  (x 

-x’)x  +  (y+y’ )y  +  (2-2*)* 
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Similarly,  without  displaying 
for  V,2G(x;x’)  is  given  by: 


V'2[G(x;x»  )a]  =  £ 


32G  {  x;  x  '  )  32G(x;x»)  32G(x;x») 

-  +  -  +  - 

3x  ’  2  3y  *  2  dz  *  2 


r325(x;x«)  32$  (x;  x  ’  )  32$(x;x')-. 

L  3x»2  3v»2  a2»2  J 


r  325  (x;  x '  )  325(x;x')  32$(x;x‘) 

L  3x*2  3 v ' 2  3  z ’ 2  J 


(31  ) 


since  a  is  a  constant  unit  vector. 


Now, 


32$  (x;  x  ’  )  r  3r  -.2  f325(x;x‘)T  30(x;x*)  32r 


3x  ’ 2 


■[£][ 


3r  3x 


(32) 


Equations  for  the  second  derivatives  of  §(x;x*)  with  respect  to 
y*  and  z’  and  the  second  derivatives  of  $(x;x’)  with  respect  to 
x’,  y’  and  z’  are  all  similar  to  Eq  (32). 

The  component  parts  of  Eq  (32)  are 


(x-x'  ) 


(33  ) 


1 


(x-x*  )2 


and 


*$<x;  x*  ) 


=  [  ik  -  -  ]  e  (x;  x  '  ) 


(35) 


and,  finally, 


a2$(x;  x*  )  r  1  -i2  5(x;  x»  ) 

-  =  I  ik - 5  (x;  x  ’  )  +  - - - 

3r*2  L  p  J  p2 


(36) 


Putting  Eqs  (33)  through  (36)  into  Eq  (32)  and  simplifying  the 
result  gives: 

325(x;x’)  j-  (x-x*)2  (x-x’)2 

-  =  5  (x;  x  ’  )  -k2  -  -  3  ik  - 

3x  ’ 2  L  r2  n3 


(x-x*  )2  i  ik  -| 

+  3  - - -  - 

p4  r2  p  J 


(37  ) 


S  imi  1  arl  y, 

32} (x;  x ’  )  r  (y-y')2 

-  =  0  (x;  x  ’  )  -k2  - 

3y’2  *-  r2 

(y-y*  )  2 


(y-y’  )2 


1  ik  -| 

*  7~  J 


(36) 
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and 


32$(x; x»  )  p  (z-z '  )2  (z-z  '  )2 

- - -  =  $  (x;  x  ’  )  -k2  -  -  3ik  - 

*  2  L  p2  p 3 


( z-z  ’  )' 


(39) 


for  the  y*  and  z*  derivatives  of  $(x;x*). 

The  second  derivatives  of  0(x;x’)  with  respect  to  x’,  y’, 

and  z*  are  nearly  identical  to  Eqs  (37  ),  (38  ),  and  (39) 

respectively.  These  derivatives  are: 

*2$(x;x’)  r  (x-x’)2  (x-x’)2 

-  =  $  (x;  x  ’  )  -k2  -  -  3  ik  - 

Av '  2  L  r2  ?3 


(x-x  »  )2 


i  i  k  -j 

-  -  +  - 

pC  p  J 


(40) 


a20(x;x’)  „  r  „  (y+y*  )2 

-  =  5 (x; x ’  )  -k2  - 

a v • 2  L  p2 


-  3  ik 


(y+y'  )2 


(y+y  )2 


i  ik 

-  -  ♦  - 

~  p  ~ 

rc  r  J 


( 4i  ) 


a2$(x;  x* ) 


5  (x;  x  *  )  £-k2 


(z-z»  )- 


(z-z  '  )' 


-  3 ik 


( z-z  ’  )2 


1  i  k  -j 

—  +  - 

rc  r  J 


(42) 
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Using  Eqs  (37)  through  (42)  in  Eq  (31)  and  noting  that 
(x-x')2  ♦  (y-y’ )2  ♦  (2-2’ )2  =  r2  and  (x-x*)2  +  (y+y* )2  ♦ 
(2-2*)  =  r2,  gives  the  final,  simple,  result  for  V’2G(x;x’) 
as: 

V* 2G (x;  x ’  )  =  -k2$  (x;  x  ’  )  -  k2$(x;x’)  =  -k2G(x;x’) 

for  x  /  x’  and  x  £  x'  (43 ) 

Eq  (43)  is  an  important  result.  It  allows  the  use  of  the 
half  space  Green's  function,  G(x;x'),  in  the  same  equations 
used  by  J.A.  Stratton  in  his  derivation  of  the  electric  and 
magnetic  field  integral  repr®? entat ions  (20:465).  The  Stratton 
equations  were  developed  with  the  free  space  Green’s  function, 

0 ( x; x '  )  which  satisfies  V’  x  V’  x  $a  =  k25»  +  V’(a  V’O) 

(20:465).  Stratton’s  results  are  based  on  this  relation  which 

is  now  also  satisfied  by  G  as  seen  by  combining  Eq  (43)  with  Eq 
(18): 

V’x  V '  x  G  (x;  x  ’  )  a  =  £k2G(x;x’)  +  V’  [a-V’G(x;  x’  )] 
for  x  X  x’  and  x  X  x’  (44) 

Using  the  results  given  in  Eq  (44)  and  the  development  of 
Eq  (3),  the  two  components  for  the  right  hand  side  of  Eq  (16) 
can  now  be  cast  in  a  more  useful  form.  In  particular, 
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G (x;  x *  )a  *  V‘x  V*  x  E(x '  )  -  E (x '  )  •  V*  x  V‘x  G (x; x *  )a 
r  a‘£iu>pJ(x’)  -  V ' x  I(x')  ♦  k2E(x’)  j 


A 

-  a 


[  K2E(X 


)  -  VG  (x;  x  *  )  [V*  •  E (x  '  )] 


] 


+  V'*  E (x '  )  £a-V*G(x;  x*  )]  (45) 

Eq  (3)  and  Eq  (45)  are  used  in  Eq  (16)  along  with  the 
following  identity: 


[V’x  I  (x’  )]  G(x;x')  dv' 


■  !hx-' 


)  x  I  (x  ’  ))  G  (x;  x  *  )  ds  ’ 


ft  A 

J(x 

t 


'  )  x  V’G(x;  x5 )  dv’ 


(46) 


An  application  of  the  divergence  theorem, 


V  •  F  dv  = 


F  .  n(x* )  ds 


(47  ) 


to  that  result  yields: 


JIf  H<s' 


)G  (x;  x  ’  )  -  (7'x  I  (x  ’  ) )  G  (x;  x  ’  ) 


i-  G  (x;  x  ’  ) 


P(x’  ) 


dv 1 


JJ[ 


iwp  [n  (x  ’  )  x  H(x '  ) )  G  (x;  x  '  )+  [n  (x '  )  x  E(x»)J  x  V»G(x;x') 


♦  [n (x ’ ) • E (x '  )J  V»G(x 


)  ] 


ds 


(48) 
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where  the  common  factor  of  the  arbitrary  constant  vector,  a  has 
been  removed.  (See  Eqs  (47),  (45),  (19),  and  (21).) 

The  integrand  of  Eq  (48)  is  singular  at  r  =  0,  (i.e. 

when  (x,  y,  i)  =  (x',y',z')).  This  violates  the  conditions  for 
use  of  the  vector  Green's  identity,  Eq  (10),  unless  the 
singular  point  is  excluded  from  the  volume,  V.  This 
singularity  actually  helps  in  the  development  of  solvable 
equations.  A  small  sphere  of  radius  rj  is  circumscribed  about 


Figure  2.  Exclusion  of  the  Singularity  by  Enclosing  It  Within 
a  Small  Sphere  of  Radius  rj  . 
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the  point  (x, y,  z  )  (see  Figure  2)  to  exclude  this  point  from 
V.  The  outward  pointing  normal  at  the  surface  of  this  sphere 
points  toward  the  center  of  the  sphere  (out  of  the  volume)  and 
along  the  radius,  r*.  For  points,  x*  on  this  small  spherical 
surf  ace 

n (x ’  )  =  r/r 

An  element  of  surface  area  on  the  sphere  is 
ds’  =  r^  s  in ( 0  )d0d$ 

with  spherical  coordinate  system  angles  0  and  $.  Using  Eq  (30) 
for  V ’ G  <  x ;  x '  )  and  taking  the  limit  as  rj  shrinks  to  zero,  the 
integral  over  the  surface  of  the  sphere  is: 

’2  n  pn  p 

[n(x»  )  x  E  (x  *  )]  x  n  (x  *  ) 

Jo  Jo  L 

K1  -|  exp  ( ikr  ) 

-  -  ik  - r  sin(0)d0d$ 

rl  -*  rt  1 

♦  0(ri )  (49) 

which  yields  4n  E(x)  for  vanishing  Rt,  since 

E(x’)  =  [n(x’)  x  E(x’)J  x  n(x’)  +  [n  (x  '  )  •  E  (x  ’  )  ]  n  (x  •  ) 

Thus,  the  contribution  of  the  volume  integral  on  the  right  hand 
side  of  Eq  (48)  is  4x  E(x)  due  to  the  volume  integral  about 
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the  singular  point.  The  electric  field  at  the  point 
x  =  (x, y, 2)  in  the  half  space  above  the  scatterer  is  given  by: 


1  rrr  r 

E(x)  =  .  .  iwpJ  (x  ’  )G  (x;  x 

4it  JJJ  L 


*  )  -  I (x '  )  X  V’G (x;  x* ) 


P(x*  ) 


+  V’G(x;  x*  ) 


] 


dv ’ 


iu,p[n(x’  )  x  H  (x  ’  )  ]  G  (x;  x 1  ) 


+  [n(x*  )xE(x' )]xV'G(x;  x' )  +  [n  (x ’  )  •  E(x ’  )  JV»G  (x;  x '  )J  ds  ' 

(50) 

Note  that  due  to  the  singularity,  the  integral  over  the 
volume  must  be  evaluated  in  the  principle  value  sense,  meaning 
that  a  vanishingly  small  spherical  volume  centered  at  x=^’  must 
be  excluded. 

With  the  choice  of  the  half  space  Green’s  function  as  the 
kernel  in  Eq  (16),  the  resulting  equation  for  E(x)  when  x 
approaches  the  surface  of  the  volume  is  the  same  as  Eq  (50), 
provided  the  surface  integral  is,  too,  a  principal  value 
integral.  The  volume  integral  remains  unchanged  since  a  small 
sphere  can  always  be  c ircumscr ibed  about  the  point  x  as  it 
approaches  S.  For  the  surface  integration  in  Eq  (50),  a  result 
due  to  Yaghjian  (23,24)  is  employed: 
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1  im 
x  — S 


V’$(x; x’  )ds ’  +  2wn (x ’  ) 


(50a  ) 


Z'  f1 


S 


V*  5  (x;  x  *  )ds  ’ 


Z'Z' 


S 


where  the  +  (-)  sign  corresponds  to  approaching  S  from  the  side 
into  (out  of)  which  n(x)  points  and  the  principal  value 
integral  employs  a  vanishingly  small  circular  patch.  For  the 
half  space  Green’s  function,  G(x;x’),  contributions  with  both 
signs  cancel,  leaving  Eq  (50)  intact  for  x  approaching  S,  with 
all  integrals  evaluated  in  the  principal  value  sense. 

A  similar  procedure  may  be  used  to  find  an  equation  for 
the  magnetic  field,  H(x),  where  H(x’)  is  used  for  the  P  vector 
of  Eq  (10);  an  equation  similar  to  Eq  (50)  may  be  developed.  A 
decidedly  simpler  procedure  is  to  use  the  concept  of  duality 
with  Eq  (50)  (4).  Based  on  the  symmetry  of  Maxwell’s 

Equations,  Table  I  shows  the  quantities  to  be  interchanged  from 
Eq  (50)  to  get  an  equation  for  H(x). 


Table  I.  Duality  Variable  Interchanges 


Quantities 

From 

To 

E(x) 

H  ( x  ) 

E 

H 

H 

-E 

P 

e 

P 

-  P 

m 

J 

-  I 

I 

-  J 

(4:  98,  99  ) 
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Using  Table  1,  the  H(x )  equation  becomes: 


H(x)  = 


-bffil 


iwel  (x  ’  )G  (x;  x  ’  )  +  J(x’)  x  V,G(x;x‘) 


P  (x>  ) 


+  V’  G  (x;  x  ’  ) 


]  dv. 


•  )  x  E  (x  ‘  )]G(x;  x*  ) 


-  [n(x’ )xH(x*  )JxV*G(x;  x*  )  -  [n (x *  )  •  H(x *  ) ]7 ’ G (x;  x ’  )j  ds  ’  (51) 

Eqs  (50)  and  (51)  may  be  further  simplified  in  certain 
problems  by  noting  that  the  volume  integrals  yield  an  incident 
field  for  a  scattering  problem  and  by  defining  surface  current 
densities,  K(x')  and  H(x’)i  and  surface  charge  densities,  a(x’) 
and  om(x'  )',  based  on  the  outward  pointing  normal  for  each 


region  considered: 


I(x’  )  =  n (x ’  )  x  H ( x 1  ) 


M ( x  ’  )  =  n  (x  ’  )  x  E  (x  1  ) 


CTm(x 1 )/p  =  n (x ’ )  •  H ( x ’ ) 


a (x ’  )/e  =  n (x ’  )  •  E  (x  ’  ) 


Using  the  simplifications  noted  above  leaves: 


(52  ) 


(53  ) 


(54) 


(55) 


E  (x )  = 


bill 


i  wyK(x  ’  )G  (x;  x  ’  )  +  H(x’)  x  V*G(x;x') 


<T(X’  ) 


V’G (x; x’  )  ds  »  +  Einc (x ) 


x  not  on  S 


(56) 
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and 


1 

H(x)=  — 
4it 


x  not  on  S 


iu>cM(x  *  )  G (x;  x  1  )  -  K(x’)  x  V‘G(x;x») 


o  (x  *  ) 

-  V  *  G (x;  x • )  da  ’  +  Hinc(x) 

p  J 


(57  ) 


where  the  'superscript  "inc"  indicates  and  incident  field. 

Eqs  (56)  and  (57)  are  valid  for  all  points  within  the  volume  in 
the  upper  have  space  of  Figure  i.  When  the  integrals  are 
interpreted  in  the  principal  value  sense,  these  equations  are 
also  valid  for  points  on  the  bounding  surface.  However,  this 
surface  includes  not  only  the  plane  between  the  upper  and  lower 
half  spaces,  but  also  includes  the  semicircular  "cover",  SR,  of 
the  upper  half  space.  Fortunately,  the  Sommerfeld  radiation 
conditions  (30)  make  the  integrals  over  SR  vanish.  The  results 
for  x  approaching  S  are: 


E(x)  = 


-i  PP  r 

JJ  L 


iupjj  (x 1  )G  (x;  x  ’  )  +  M(x  *  )  x  V’G(x;x') 


S 


<MX  ’  ) 

♦  -  V 


’  G(x;  x  *  )  J 


ds  ’  +  Einc  (x  ),  x — *S 

(58) 


and 


H(x)  = 


Hi 


iweM(x’)  G  (x;  x  ’  )  -  K(x’)  x  y'G(x;x’) 


°  (S’  )  -) 

* -  V '  G  (x;  x  ’  )  ds  ’  +  Hinc(x),  x-»S 

p  J 


(59) 
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Eqs  (56)  through  (59)  are  the  equations  which  represent 


the  total  electric  or  magnetic  field  at  any  point  inside  or  on 
the  surface  of  the  upper  half  space  of  Figure  i  in  terms  of  the 
surface  charges,  ct  and  om,  and  the  surface  currents,  K(x’)  and 
M(x’).  To  find  these  currents  and  charges,  equations  must  be 
developed  to  find  the  fields,  currents,  and  charges  inside  the 
cavity  of  Figure  1  because  the  two  spaces  are  coupled  together 
at  the  interface  between  the  surface.  Recall,  this  interface 
has  been  named  the  aperture. 

Ins ide  the  Cavity 

Inside  the  cavity,  the  free  space  Green’s  function  in  the 
a  direction  may  be  used  for  Q  (20)  and  E  or  K  used  for  P  in  Eq 
(10).  The  volume  in  the  integral  is  the  volumi  of  the  cavity 
and  the  surface  is  the  surface  of  the  cavity.  Again,  no 
sources  are  found  within  the  cavity  and  the  surface  currents 
and  charges  may  be  modeled  as  they  were  in  Eqs  (52)  through 
(55).  The  fields  incident  to  the  aperture  will  induce  surface 
currents  and  charges  across  the  aperture.  These  surface 
currents  and  charges  will  generate  the  fields  that  are  incident 
to  the  cavity.  In  this  way,  the  cavity  is  coupled  to  the  half 
space  above  the  aperture.  Using  all  of  the  above  information, 
the  electric  and  magnetic  fields  inside  the  cavity  are 


E(x)  = 


-i  rr  r 

4tt  JJ  L 


iu»pK(x  ’  )$  (x;  x  ’  )  ♦  M(x’)  x  V’$(x;x*) 


a(x’  ) 


V'$(x 


x;x'  )  ] 


ds 


for  x  not  on  the  surface,  S 
and 


(60) 


H (x  )  - 


4ir 


iweM(x’)  $(x;x’)  -  K(x’)  x  V’$(x;x') 


cr  (x  *  ) 

-i&-= —  V'  5(x;  x'  )  ds  » 

IJ  J 


for  x  not  on  the  surface,  S 


(61  ) 


To  let  the  observation  point  approach  the  surface  S  bounding 
the  cavity,  Eq  (50a)  is  again  employed.  The  results  are: 


-1 


E(x) 


4w 


ff( 


iwpK  (x  *  )  0  (x;  x  ’  )  +  M(x’)  x  V’0(x;x’) 


a(x’  )  -i 

+  -  V’5(x;x’)  ds  ’ 

e. 


1  i  O(X) 

—  M(x )  x  n(x)  - - n  (x  ’  ) 

2  2  e 


for  x  -*  S 


(62a) 
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and 


HJ[ 


H (x  )  =  -  -H-  |  iu»eH(x  ’  )  5  (x;  x  ’  )  -  K(x'  )  x  V'$(x;  x*  ) 

4tt 

S 


a  (x’  )  1 

—  V’  5  (x;  x  *  )  ds 
U  J 


1  i  CTm(X) 

~  -  ?  (x  )  x  n  (x  )  - - n  (x  ’  ) 

2  2  e 


for  x— S 


(62b) 


With  Eqs  (56)  through  (62)  the  fields  at  any  point  x  can 
be  determined  if  the  incident  fields  and  the  surface  charges 
and  currents  are  known.  The  next  step  in  the  derivation  of  the 
integral  equations,  therefore,  deals  with  finding  the  surface 
currents  and  charges. 

Find ing  the  Currents  and  Charges 

The  currents  and  charges  needed  to  use  Eqs  (56)  through 
(62)  are  found  by  the  solution  of  a  coupled  set  of  second  kind 
Fredholm  integral  equations  which  will  be  derived  presently. 

To  help  in  the  computations,  Eq  (62a)  is  multiplied  by  iu»e . 

This  allows  for  the  solution  of  iweH(x’)  instead  of  solving 

for  M(x’)  alone.  Notice  that  iu>e  multiplied  by  iwp 

yields  -k2  by  Eq  (8).  It  is  also  more  convenient  to  let  the 
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constant  1/4*  be  part  of  the  Green's  functions  so  that 
$(x;x')/4w  is  referred  to  everywhere  as  ${x;x').  Also,  let 
ffjn/Mi  from  here  on,  be  referred  to  as  om. 

In  the  equations  to  follow,  Sc  refers  to  the  surface  of 
the  cavity  excluding  the  aperture  and  A  refers  to  the  aperture; 
as  before,  S  refers  to  the  surface  of  the  upper  half  plane 
(which  also  includes  the  aperture).  For  the  magnetic  current 
across  the  aperture,  Eqs  (53)  and  (62a)  along  with  the  above 
simplifications  yields  a  second  Rind  Fredholm  integral  equation 
for  the  current  H  forced  by  integrals  of  the,  as  yet,  unknown 
current  K  and  charge  a: 


1 

—  iu>€M(x) 
2 


a  .  « 

n  (x )  x 


E(x> )$ (x; x ’  ) 


-  M(x’  )  x  V’$(x;x’)iwe 


where  x  €  A. 


-  i<DCT  (x  ’  ) 


V’ $ (x;  x  ’  ) 


(63  ) 


Instead  of  using  Eq  (62a),  Eq  (56)  could  be  used  to  find 
Eq  (63).  However,  since  Eq  (58)  is  for  the  upper  half  space, 
the  half  space  Green's  function  would  be  used  in  place  of 
$(x;x')  and  the  tangential  component  of  the  effective  incident 
electric  field  across  the  aperture  would  be  added.  Eq  (62a)  is 
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chosen  in  place  of  Eq  (58)  since  the  former  involves  an 
integration  over  the  bounded  aperture,  whereas  the  latter 
requires  integration  over  the  unbounded  xz  plane. 

A  further  simplification  is  possible  since  M(x)  =  0  when  x 
is  on  a  perfect  electric  conductor  (4:34).  The  cavity  walls 
are  assumed  to  be  PEC.  Furthermore,  to  avoid  confusion 
regarding  the  normal  used  in  defining  quantities  in  Eqs  (52) 
through  (55),  the  magnetic  current  is  defined  now  to  be 
M(x’)  =  -y  x  E(x’).  where  y  is  the  unit  vector  pointed  in  the 
positive  direction  along  the  y-axis.  Using  this  definition, 
the  equation  for  the  magnetic  current  when  x  is  on  the  aperture 
becomes: 


1  PP  f  1 

—  iweM(x)  =  - k2  1^  y  x  |(x'  )}(x;  x'  )  J  ds  ’ 


Sr  U  A 


-  iwe  + 


I5 


y  x  H(x’)  x  V’5(x;x’)  dx  ’ 


+  It  y  x  i<i><T(x’)  V’5(x;x’)  ds’ 


Sc  U  A 


x  €  A  (64) 


To  find  the  electric  current  across  the  aperture,  Eqs  (52) 
and  (59)  will  be  used  with  the  current  defined  as 
K(x’)  =  y  x  H(x’)-  The  result  is 
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I(X)  =  y  x  Hinc(x)  +  y  x  j-j-  iweM(x  ’  )G  (x;  x  *  ) 

S 

-K(x’)  x  V*  G  (x;  x '  )  -  cr  V’G(x;x‘)  1  d  s 
-  __  m  - J 


Jf[ 


with  x  €  A 


(65) 


(Recall  that  cm/y  was  earlier  redefined  as  cm  for  simplicity  of 
presentation. ) 

Although  this,  too,  requires  integration  over  the 
unbounded  xz  plane,  some  means  of  coupling  the  incident  field 
into  the  problem  must  be  employed.  It  will  be  explained  in  the 
next  chapter  how  the  judicious  choice  of  Green’s  function, 
G(x;x’),  combined  with  other  problem  specifics  is  employed  to 
eliminate  the  integrations  over  the  plane. 

For  the  electric  current  on  the  surface  of  the  cavity 
(minus  the  aperture),  Eq  (62b)  is  used  in  conjunction  with  Eq 
(52).  The  equation  for  the  electric  current  in  the  cavity 
becomes: 


1 

-  I(x) 
2 


ff[ 


n  (x  )  x  ft  -iu»eH(x’  )}  ( x;  x’  )  ♦  I(x’)  x  V’5(x;x’) 


Sc  U  A 


or  V*  $  (x;  x  ’  )  1  di 
m  - J 


with  x  €  Sc 


(66) 
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To  find  an  equation  for  the  electric  charge  density  on  the 
surface  w 2  the  cavity  including  the  aperture,  Eq  (62a;  is  used 
in  conjunction  with  Eq  (55).  Since  a  is  always  multiplied  by 
iw  in  Eqs  (58)  through  (62),  it  is  appropriate  to  find  an 
equation  for  i wct.  Again,  M(x)  ;  0  when  x  is  on  the 
surface  of  a  PEC,  and  n (x )  is  dependent  upon  the  cavity 
surface.  The  equation  for  icja  is: 


—  itdcr(x) 
2 


JJ 

A 


a  / 

n  (x 


[m(x 


)  x  V*  5 (x 


*!*’>] 


ds 


+  4  f  n  (x  ) 


Sc  U  A 


■>  p 

A  , 

-  ■  ■  n (X  ) 

J  L 

Sc  U  4 
with  x  €  Sc  U  A 


K2  K (x ’  )  0 ( x ;  x ’  )  ds  ’ 


iwer  (x  ’  )  V’  {  (x;  x  ’  )  ds  ’ 


(67  ) 


An  analogous  equation  for  am  could  be  obtained  from  Eq 
(62b)  but  is  not  listed  here.  Eqs  (64)  through  (67)  and  the 
equation  for  <7m  are  the  necessary  equations  for  the  surface 
current  and  charge  densities.  Each  is  a  second  kind  Fredholm 
integral  equation.  The  equations  are  coupled  and  may  be  viewed 
as  a  linear  operator,  L,  acting  on  the  vector  of  unknowns, 

X  =  [M(x),  K(x),  cr(x),  am(x  )],  according  to  LX  =  B,  where 
B  =  (0,  y  x  Hinc(x),  0,  0).  The  operator,  L,  is  a  matrix  of 
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linear  operators  having  second  kind  Fredholm  operators  along 
the  d  i  2  «  ~  r.  a  1  .  Alternate,  simpler,  f  ormu !  2*  i  ons  v.'hich  exciud- 
charge  densities  are  possible.  In  one  such  formulation,  the 
resulting  operator  involves  a  first  kind  Fredholm  operation, 
inversion  of  which  is  possibly  ill-posed  (2).  Another 
formulation  replaces  iuicr  by  Vs  •  K  and  iu>cjm  by  -Vs  ■  H  in 
which  Vs  is  the  divergence  operator  in  the  surface.  These 
relations  follow  directly  from  Maxwell’s  equations  and  are 
commonly  employed  (15)  when  using  surface  integral  equations. 
However,  the  formulation  herein  is  chosen  to  avoid  numerical 
differentiation  for  approximate  solutions  obtained  in  a  later 
chapter . 
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II!  . 


Seduction  tn  a  Two-d  <mens  ion  at  Prob  1  eg 


The  Cavity 

For  this  problem,  a  two-dimensional  cavity  (Figure  i)  of 
arbitrary  shape  is  considered.  Along  the  z-axis,  the  cross- 
sectional  shape  of  the  cavity  is  constant.  For  such  a  problem, 
the  currents  and  charges,  M,  K,  a,  and  am,  as  well  as  the 
resulting  fields,  are  independent  of  z. 

This  simplification  eliminates  integration  in  the  z 
direction  from  -a>  to  +a>  for  the  surface  currents  and  charges. 

It  also  allows  the  use  of  the  half  space  Green’s  function  in 
two-dimensions  for  G(x;x’)  and  the  two-dimensional  free  space 
Green’s  function  for  $(x;x’).  In  two-dimensions  these  Green’s 
functions  are  ( 13  ): 

G  (x;  x  ’  )  =  it  iH^1  )  (Kr  )  +  *iHjl>(kr)  (68) 

and 

$  (x;  x  ’  )  =  iriH^1  >  (Kr)  (69  ) 

where,  now 


ii/2 

(x-x ’  )2  + 

(y-y’  )2 

] 

l‘/2 

(70) 

(x-x’  )2  + 

(y+y’ )2 

(71  ) 

and 

H^1)(kr)  is  a  Hankel  function  of  the  first  kind  of  order  zero. 
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Dividing  Eqs  (68)  and  (69)  by  4w  will  help  simplify  Eqs  (58) 
through  (61).  Thus,  the  Green’s  functions  are  now  redefined  to 
be 


G(x;  x’  ) 


and 


i  i 

-  Hj1 ) (Kr)  +  -  Hjl ) (Kr) 

4  4 


(72) 


$  (x;  x ’ ) 


i 

4 


H^1 > (Kr) 


(73  ) 


In  two-dimensions,  the  surface  integrals  of  Eqs  (58) 
through  (61)  are  now  line  integrals  for  finding  the  "surface" 
currents  and  charges. 

The  Inc  ldent  Fields 

Assume  that  the  field  originates  at  a  long  distance  from 
the  aperture  such  that  as  it  impinges  upon  the  aperture  the 
field  is  a  plane  wave.  Further,  let  the  incident  field  be  a  TE 
(Transverse  Electric  Field)  plane  wave.  That  is,  the  electric 
field  is  transverse  to  the  z  direction  at  all  times.  Also,  let 
the  medium  above  the  plane  be  free  space  (this  has  already  been 
assumed)  and  let  both  the  electric  and  magnetic  fields  always 
be  transverse  to  the  direction  of  propagation.  This  is  a  TEM 
(Transverse  Electromagnetic)  plane  wave.  For  a  TE  plane  wave 
traveling  in  the  TEM  mode,  the  H  field  is  in  the  z  direction. 
Referring  to  Figure  3,  0A  is  the  positive  angle  measured 
clockwise  from  the  y  axis  along  which  the  wave  propagates. 

Such  an  incident  plane  wave  with  exp(-iwt)  time  dependence 
(assumed  at  the  outset)  has  spatial  dependence  governed  by 
(4:  146  ): 
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With  unit  magnitude,  the  incident  field  is: 


jjinc(x)  .  J  exp{-ik[x  sin(6^)  +  y  cos(0i)})  (75) 


With  this  incident  field  and  the  two-dimensional  cavity 
discussed  above,  Eqs  (64)  through  (67)  may  be  simplified.  The 
simplifying  observations  are: 


1.  From  Eq  (54),  the  magnetic  surface  charge  density,  crm(x), 
is  zero  for  the  TE  plane  wave  considered  above. 

2.  On  the  PEC,  the  magnetic  surface  current  density,  H(x),  is 
zero . 

3.  The  gradient,  in  primed  coordinates,  of  the  Green’s 
function,  V’G(x;x’),  has  only  an  x  component  when  either  y  or 
y’  equals  zero.  This  is  the  Neumann  boundary  condition 
discussed  in  the  introduction  of  the  thesis.  Therefore, 

K(x’)  x  V'G(x;x’)  =  0  when  x’  is  on  the  x-axis. 

4.  The  gradient  of  the  Green’s  function  never  has  a  z 
component  in  this  two-dimensional  problem.  Therefore, 

H ( x ’ ) • V ’ G ( x;  x ’  )  =  0  for  all  x’. 

5.  The  three-dimensional  surface  is  now  a  two-dimensional 
contour.  The  equations  to  follow  will  reflect  that;  thus,  Sc 
now  indicates  the  contour  around  the  cavity  and  A  indicates  the 
line  across  the  aperture. 


Using  these  simplifications,  the  five  equations  that  will  be 
approximated  in  the  following  analysis  are: 


H (x  ) 


Hinc (x) 


A 


iweM(x ’  )G(x; x’  ) 


dx  ’ 


(76) 
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1 

—  ru>eH(x) 
2 


Sc  l 

with  x  € 


K  (x )  =  y  X 


with  x  € 


1 

-  I(x)  =  - 
2 


with  x  € 


k2 


Sc  U  A 


K (x '  )  5  (x;  x  ’  ) 


_ 


dc  ' 


y  x  iu»CT(x’)  V’0(x;x’)  dc* 
A 


A 


(77  ) 


Hinc(x)  +  y  x  -l  iwcM(x  ’  )G  (x;  x  ’  )  dx  ’ 


A 


A 


(78) 


*  /  « 

n  (x  )  x 


+ 


J-  iweH(x  ’  )0  (x;  x  ’  ) 

A 

j-  K(x»  )  x  V’  0  (x;  x  •  ) 
U  A 


dx  ’ 


dc  * 


SC 


(79) 
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and 


1 

—  lliKT(X) 
2 


-  me 


» 

<-  n  (x  )  •  ^Mfx  '  )  x  V’  $ (x; x 


+ 


n ( x ) • K2 


U  A 


K  (x  ’  )  5  ( x i  x  ’  ) 


dc  ’ 


•  •  n  (x  )  •  iu>a  (x  ’  )  V’$(x;x>)  dc  ’ 
Sc  U  A 

with  x  €  Sc  U  A 


dc  ’ 


(80) 


Eqs  (77)  through  (80)  are  a  coupled  set  of  Fredholm 
integral  equations  of  the  second  kind  for  finding  the  currents 
and  charges  needed  to  solve  for  the  scattered  fields,  using  Eq 
(76),  from  the  two-dimensional  cavity  introduced  at  the 
beginning  of  the  chapter.  Recall  if  iwcr  is  replaced  by 
Vs  • K,  as  described  at  the  end  of  Chapter  2,  Eq  (80)  is  not 
needed.  In  either  case,  the  equations  are  specifically 
designed  for  use  with  the  TE  plane  wave  of  Eq  (75).  Solutions 
of  these  equations  may  be  approximated  numerically  by  any  of 
several  methods.  The  chosen  method  of  approximating  the 
solutions  to  Eqs  (77)  through  (80)  is  introduced  in  the  next 
chapter . 
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IV.  Approximate  Solution  of  the  Integral  Equations 

In  Eqs  (77)  through  (80),  linear  operations  (the  integrals 
and  vector  operations)  act  on  a  set  of  unknowns  (the  currents 
and  charges)  and  must  be  equal  to  a  set  of  Rnovns  (the  incident 
fields).  If  the  inverse  of  the  linear  operators  can  be  found 
and  applied  to  the  set  of  knowns,  the  set  of  unknowns  can  be 
determined.  The  Method  of  Moments  is  a  reasonably  simple  means 
by  which  one  can  approximate  such  an  operator  equation  for 
numerical  solution  (3).  Eqs  (77)  through  (80)  can  be  solved 
using  moment  methods  by  using  the  following  approximations: 

M 

iweM(x)  ~  E  an  f n (x ) 
n  =  1 

N 

E(x)  2  E  Pn  hn(x) 
n  =  1 

N 

iujq(x  )  3  E  rn  pn  (x  )  (8!  ) 

n  =  i 

where  an,  0n,  and  rn  are  constants  and  fn(x),  Pn(*)>  and  hn  ) 
are  functions  to  be  described  below. 

The  "3"  sign  in  Eq  (81)  implies  approximations  because  the 
numbers  M  and  N  are  assumed  to  be  finite.  As  M  and  N  approach 
infinity,  the  summation  (right  hand  side)  approaches  a  smooth 
curve  (left  hand  side)  and  the  "2"  can  be  replaced  by  an  equal 
sign.  In  Eq  (81),  fn(x),  hn(*)*  and  pn(x)  are  basis  functions 
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of  position  (3:42).  The  functions  fn(x)  and  hn (x )  are  vectors 
because  M(x)  and  K(x)  are  vectors.  Using  Eq  (81)  and  the 
linear  operators  in  Eqs  (77)  through  (80)  yields: 


1  1  H  H  N  N 

-  iu)eM(x)  x  -  Eanf  n  (x )  =  EanLif^x’)  +  E0nL2hn(x')  +  ErnL3pn(x') 

2  2n=l  n=l  n=l  n=l 


with  x  €  A  (82) 


N  M 

K(x)  z  E  0n  hn(x)  =  EanL4fn(x')  +  Hinc(x) 
n= 1  n=  1 

with  X  €  A  (83  ) 


1  1  N  M  N 

?  ( x  )  ~  E  Pn  ]2n  )  ~  EajjLgfj^x  '  )  +  ^*Pn^*6l?n  '  ) 

2  2  n  = 1  n  =  l  n=l 


with  x  €  Sc  (84) 


1  1  N  M  N  N 

-  iwa(x)  z  -  Ernpn(x)  =  EanL7fn(x'>  -  E0nL8hn(x')  +  ErnL9pn(x’) 

2  2  n  =  1  n  =  1  n=l  n:l 

with  x  €  Sc  U  A  ( 85 ) 

The  linear  operators  through  L9  are  all  functions  of 
position.  That  is,  =  Lj(x),  L2  =  L2(x),  etc.  Further,  the 

linear  operators  may  be  determined  by  examining  the  equations 
from  which  they  came:  Eqs  (77)  through  (80).  For  instance, 
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L  f  (X’  ) 
i_n  _ 


r< 


y  x  f  (x’  )  x  7'{(x;x'  )  dx’ 

_  n  _  _  _ 

A 

The  method  is  employed  by  defining  a  weighting  function, 
Wm(x;xm).  for  the  specific  problem  being  solved,  and  taking  the 
inner  product  of  Wm  with  each  of  Eqs  (77)  through  (80);  where 
the  inner  product  is  defined  as: 


/  (x,  x  ),  f  (x  )>  =  - -W  (x,  x  )f  (x)  dc  ' 
m  ~  “m  mm 

i 


(86) 


This  yields 


MM  N 

1/2  < Wjjj,  f  n  ( x  )  >  =  Ean<wm*  iJ n  ^  +  ^Pn  ^^m>  ^2— n  ( JS  *  )  ^ 

n=l  n  =  1  n=l 


N 

+  ETn<Wm,  L3pn(x*  )> 
n  =  1 

with  x  €  A,  ( 87  ) 


N  M 

EPn<wm«!]n<2S>>  =  Ean<Wm,  L4f  n  (x  ’  )  >  ♦  <Wm,  Hinc(x)> 
n  =  1  n  =  1 


with  x  €  A, 


(88  ) 
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and 


N  M  N 

1/2  E0n<Wm,  hn  (x  )  >  =  Ean<Wm,  L5J n  (x  1  )  >  ♦  E0n<Wm,  Lghn  (x  ’  )  > 
n=l  n=  1  n=l 

with  x  €  Sc,  (89) 

N  M  N 

1/2  E  y  <  W|jj,  Pn  ( Xy  )  >  -  £an<Wm,  L7^n(x’)^  +  ^Pn  ^^m»  (£ 1  )  ^ 

n  =  1  n  =  1  n=  1 

N 

+  Ern<Vm,  L9pn  (x  ‘  )> 

n  =  1 

with  x  €  Sc  U  A.  (90) 

Eqs  (87)  through  (90)  allow  the  currents  and  charges  to  be 
found  using  moment  methods  directly.  No  matter  what  the  shape 
of  the  cavity’s  cross  section,  Eqs  (87)  through  (90)  should 
yield  a  reasonably  accurate  solution  if  M  and  N  are  chosen 
judiciously.  However  it  is  advantageous  to  deal  with  a 
specific  cross  section  as  opposed  to  a  generic  or  arbitrary 
cross  section.  For  this  reason,  the  specific  cavity  cross 
section  of  interest  is  introduced  in  the  next  section. 
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On  Figure  4,  the  outward  pointing  normals  are  indicated  on 
each  surface.  Further,  the  electric  current  direction  is  also 
noted.  The  integrations  within  the  linear  operators  of  Eqs 
(87)  through  (90)  are  line  integrals  with  respect  to  positive 
arc  length  around  the  cavity  and  across  the  aperture  (21). 

Table  II  gives  a  list  of  variables  normalized  to  the 
width.  The  first  column  gives  the  original  variable.  The 
second  column  gives  the  normalized  variable  and  the  third 
column  gives  the  name  used  in  the  computations. 


Table  II.  Quantities  Normalized  to  the  Width 


|  Original 

Variable 

Normal ized 
Variable 

Used 

Quantity 

w 

w/w 

i 

X 

x/w 

X 

y 

y/w 

y 

z 

z/w 

z 

d 

d/w 

d 

X 

X/w 

X 

k 

kw 

k 

With  the  geometry 
broken  up  into  N  equal  1 


indicated  in  Figure  4, 
y  sized  segments  with 


the 
M  of 


cavity  can  be 
them  across 
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the  aperture.  If  the  size  of  each  segment  is  termed  ACn 
(actually  normalized  to  w)  where  AC}  =  AC 2  =  ...  = 

AC)],  the  basis  functions  can  be  defined  as: 

r  le  if  x  is  on  AC 

t  (x)  =  "  n  (91) 

n  L  0  elsewhere 

where  the  direction,  e,  associated  with  fn(x)  is  the  same  as 
that  of  the  current  for  each  particular  segment.  For  instance, 
on  the  inside  of  the  aperture  J(x)  is  in  the  k  direction. 
Similar  basis  functions  are  defined  for  hn  and  pn. 

Using  "point  matching"  for  testing,  the  weighting 
function,  Wm(x;x'),  is  defined  as  (3:42) 

Wm(x;x»)  =  W(x;xm)  =  $(x  -  xm)  (92) 

Eq  (91)  reduces  the  integrals  within  the  linear  operations 

L| (x  ),  L2 (x  ) . L9 (x )  to  integrations  over  a  small  segment, 

AC.  Eq  (9*)  replaces  the  position  vector  x  with  the 
position  vector  xm  in  Eqs  (87)  through  (90).  To  illustrate, 
let  r(x;x’)  be  any  function  that  is  continuous  along  the  line  L 
and  let  the  line,  L,  be  broken  into  H  segments.  Then 


> 

K 


r (x;  x’  ) 


dx  * 


M  p 

E  r (x; x’  ) 
n  =  1  v 


ACn 


dx  • 
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Taking  the  inner  product  as  defined  above  yields 


f  [ J 


r  (x;  x  ’  )dx 


S(x-x  ) 

m 


dx 


AC 


n 


M  r 

E  r(x  ; x* )  dx* 

n=  i  J  “ 


If  r(xjx')  is  one  of  f,  h,  or  p  given  in  Eq  (9i),  only  those 
line  segments  ACn  survive  on  which  r  is  non-zero. 

Referring  now  to  the  problem  at  hand,  by  letting  xm  be 
the  midpoint  of  each  line  segment,  ACm,  (as  m  goes  from  1 
to  N)  and  xn  be  the  midpoint  of  each  line  segment  ACn  (as  n 
goes  from  1  to  H),  the  integrals  over  the  line  segment  ACn 
may  be  approximated  by: 

f(x  ; x  *  )  dx  *  z  f(x  ; x  )  AC  (94) 

J  _m  _  _  _m  _n  n 

ACn 

if  ACn  is  sufficiently  small. 

Using  Eqs  (91)  through  (94)  and  performing  the  inner 
products  of  Wm  with  the  linear  operators  operating  upon  their 
appropriate  functions,  Eqs  (87)  through  (90)  may  be  cast  in  a 
more  easily  solvable  form. 

At  either  y  or  y'  equal  to  zero,  G(xm;xn)  s  2$(xm;xn). 
This  fact  is  used  to  simplify  the  equations.  It  is  important 
to  note  that  the  integrations  are  line  integrals  and  are 
independent  of  direction  across  the  small  segment,  AC. 
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Since  ACn  is  always  a  positive  number,  the  line  integration 
does  not  affect  the  sign  of  the  approximation  as  defined  in  Eq 
(94).  Further,  the  principle  value  of  the  integral  over 
V’O(Xmi^n)  is  identically  zero  (as  will  be  shown)  when  n=m. 
When  n=m,  the  segment  ACn  is  termed  a  self-patch.  The 
principle  value  integrals  over  the  self-patches  will  now  be 
examined  in  more  detail. 

Self-Patch  Integrations 

In  two-dimens  ions,  V'jHxjx1)  is 

V’5(x;  x’  )  =  —  x  +  —  y 
3x  ’  dy' 

-  H(&>(kr)  r  "I 

=  — 1 -  j^(x-x')*  +  (y-y')£  J  (95) 

In  Eq  (95)  the  Hankel  function  of  the  first  kind  order  one  and 
the  i/r  term  both  produce  s ingularit ies  when  x  =  x*.  Performing 
the  integrals  in  the  principle  value  sense,  however,  does  not 
produce  singularities. 

When  integrating  V'5(x;x,),x  along  the  y-axis,  the  term 
involving  x-x'  (or  xm-xn  in  the  approximation)  will  always 
multiply  the  integral  by  zero  when  x=x’  (or  xm=xn).  The  same 
situation  is  involved  when  integrating  V’5(x;x’)'y  along  the 
x-axis,  where  the  (y-y* )  term  always  multiplies  the  integral 
by  zero  when  y  =  y*.  Thus,  those  particular  self-patches  are 
identically  zero. 
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For  the  self-patches  involving  integrations  in  the  x 
direction  of  V'$(x;x*)'X,  a  more  rigorous  proof  is  necessary. 
The  same  situation  applies  when  integrating  in  the  y  direction 
of  V’{(x;x'  )7-  In  the  following  proof,  the  variable  "x"  may 
be  replace  by  the  variable  My"  to  achieve  the  same  result  as 
long  as  every  "x"  is  replaced  by  a  "y"  . 

When  integrating  on  a  self-patch  in  the  x  direction  (for 
the  rectangular  channel),  y  is  always  equal  to  y*  .  Thus, 


xm+AC/2 


4  V  *  $ (x  ;  x  *  ) • x  dx 

m 

*j 

ACm 


i  H<‘ ) [k I x  -x'  I) 

*  :  —  --  - 1 - m -  (X  -X  ’ 


IV*' I 


m 


)  ax  * 


xm-AC/2 


(96) 


Eq  (96)  is  a  symmetrical  integration  of  an  odd  function  about 
the  singularity  x*=0.  This  integration  is  therefore 
identically  zero. 

The  integrations  over  all  self  patches  involving  V’0(x;x’) 
have  been  shown  to  be  equal  to  zero.  All  of  the  self-patches 
that  do  not  involve  a  V,5(x;x*)  term  have,  instead,  only  a 
5(x;x*)  term.  These  integrals  do  not  reduce  to  zero.  Again, 
the  variable  "x"  may  be  replaced  by  the  variable  "y"  in  the 
following  integrals. 
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j  $(x  ; 

J  ■ 


xm+AC/2 

i 

x»)  dx  *  =  —  --  H< 1 ) (K I x  -x*|)  dx 
4  o  1  a  1 

Xb-AC/2 


AC/2 


:  lim  -  H*1  >(kx)  dx 

€ — mO  2  0 


(97) 


For  small  arguments,  the  asymptotic  approximation  for 
H0(1>  is  (1:360) 


H0  < 1  )  ( i ) 


Thus, 


f  S,V 


x  ’  )  dx 


1 n ( z/2  )  ♦  1 

+  i  — 

y  -  i  —  In (z/2 ) 

w 

2s 

2 

T  l2 

i  -  (l-y)  - 

1  “ 

-  +  0[2^1n(z)J 

x 

j  4 

(98) 

’s  constant 

a  0.5772156649015. . . 

KA(J  r  i 

2y  -1 

|  AC  r  -1 

S  -  - 

— 

-  —  In (kAC/4)  -  l 

4  L  k 

kw  - 

'  2s  L  -» 

r  kAC  -.3 

r  2  i 

i  r  i 

♦  - 

j  —  I 

In (kAC/4 )  -  1/3 

L  4  J 

L  3kx  -1 

1  L  J 

(99) 
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All  of  the  self-patch  integrations  have  now  been 
determined.  Those  line  integrals  over  sections  other  than 
self-patches  will  be  approximated  using  Eq  (89). 


Solving  for  the  Currents  and  Charges 

Most  of  the  equations  used  thus  far  are  vector  valued. 
However,  these  equations  can  be  broken  into  x  and  y  components 
if  necessary  and  a  set  of  equations  developed  from  each  of  Eqs 
(88)  through  (9i)  and  their  corresponding  earlier  equations, 
Eqs  (77)  through  (80). 

The  set  of  equations  corresponding  to  Eqs  (77)  and  (88) 

are: 


For  m  =  1,2,  3...  M  where  r\(xm)  =  -y 


M  M 

0  =  E  an  [V ’  0  ( XjhJ  xn  )  •  y]  ACn  -  E  Pjj  Q  (xm;  xn  )AC^ 
n=  1  n  =  1 


L+  2M 

M+L 

E  fin 

5(xm; 

Xn)ACf)  ♦  E  Tjj  [7*  0  (xm;  xn  ) 

•  X) ACn 

n  =  L+M 

n=  1 

N 

xn)-x]ACn  ♦  am/2 

E  r„  [7* 

5  (Xm> 

(100) 

n=M+L 


where 

L  =  the  number  of  segments  along  S-  and  S+ 
M  =  the  number  of  segments  along  A  and  Sd 


Eq  (100)  is  actually  a  set  of  M  equations  because  it  applies 
for  m=l,2,  ...,M.  There  are  only  M  equations  because  xm  is  not 
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on  the  aperture  when  m  >  M;  Eq  (77),  and  therefore  Eq  (80),  is 
only  good  when  the  observation  point  (x  or  xm)  is  on  the 
aperture.  The  set  of  equations  that  correspond  to  Eqs  (78)  and 
(89)  are: 

For  m  =  i,  2,  3,  .  .  .M  where  n(xm)  =  y 

1  M 

-  H*nc  (xm )  =  -  E  an  8(xm;xn)  ACn  +  pm/2  (101) 

2  n  =  i 

The  set  of  equations  that  correspond  to  Eqs  (79)  and  ( 90 )  are: 


For  m  = 

M+i,  M+2,  .  .  .M+L  where 

n(xm)  = 

A 

-x 

M 

M 

0  =  E  an 

0  (5m»  25n  J^^n  +  ^  0n 

[V 

8  (?b>  ?n ) 

•y)Acn  + 

0m/2 

It 

c 

n  =  1 

Lr  2M 

N 

-  E  0 

ntV'8(xm;  xn)-y]ACn 

+  E 

f>n  [V 

8 (xm;  xn ) 

•  X] ACn 

n=M+L 

n 

=  L+2M 

(102) 


For  M+L  <  m  <  2M+L  where  n(xm)  =  -y 


M  M 

0  -  ~E  an8  (Xmi  15n  )ACn  ~  ^  Pn  [^'{(Jb'  2n)’7l^n  "  0m/2 
n= '  n= 1 


*  N 

8n  5 (xm;  xn ) • x] ACn  -  E  8n  [V’ 0 (xm;  xn ) • x] ACn 

n=L+2M 

(103) 


L+M 
+  E 
n=M 
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For  2M+L  <  m  <  N  where  n(xa)  =  x 


M  M 

0  =  -  E  an0(xm;  xn)ACn  -  E  0n  [V> 0 (xm;  xn ) • y) ACn 
n  =  l  n=l 


-  em/2 


L+M  L+2M 

+  E  0n[V’5(xm;  xn)-x]ACn  +  E  0n  [V>3(xm;  xn)-y]ACn 


n=H 


n  =  L+M 


(104) 


And  finally,  the  set  of  equations  that  correspond  to  Eqs  (80) 
and  (91)  are: 

For  0  <  m  <  M  where  n(xm)  =  y 


M 


M+L 


0  =  E  [V 1  5  ( Xjn*  ?n  )  x  ]  ACn  +  E  [V1  {  (Xgj  Xjj  )  •  y )  AC  n 

n - 1  n  =M 


M+L 


H 


+  k  E  Pn3  (xm;  xn  )AC  -  K2  E  0n$  (im*  *n  )AC 
n=M  n  =  2M+L 


N 

-  E  Tn [V ' 0 (xm;  Xn  )  •  £] ACn  -  Tm/2 
n=M+L 


(105) 


For  M  <  m  <  M+L  where  n(xm)  =  -x 


M 


M 


o  =  E  an[V‘5(xm;  xn)-y)ACn  +  E  Tn [V’ $ (xm;  xn ) • x) ACn 
n  =  1  n-  1 


M  2M+L 

k2  E  0n6<2m5  *n)AC  +  K2  E  0n5(xmJXn)AC 


n  =  1 


n  =M+L 


N 

♦  E  Tn [V ’ 5 (xm;  xn ) • x] ACn  -  tm/2 
nrM+L 


(106) 
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For  M+L  <  m  <  L+2H  where  n(xm)  =  -y 

M  H+L 

0  =  -E  an[V’5<Xm;  xn)-^3ACn  ♦  E  rn  [V’ 8  (xm;  xn  )  •  y)  ACn 
n=l  n=l 

M+L  N 

"  k2  E  en8(xm;xn)AC  +  k2  E  3n8(xm;xn)AC 
n=H  n  =  2M+L 

H 

♦  E  rnCV'8(xm;  xn)-fr]ACn  -  rn/2  (107) 

n  =  2H+  L 


For  L+2M  <  m  <  N  where  n(xm)  =  x 

M  M+L 

0  =  ~E  On  [V '  8  (?in*  2?n  ) '  Tl  ~  E  Tjj  [V  8  (X|j!  *n ) '  *)  ACn 
n=  1  n  =  1 


M  2M+L 

+  k2  E  Pn0(XBi;  xn)AC  -  K2  E  PnC(xm;xn)AC 
n  = 1  n  =M+  L 


2M+L 

-  E  rn [V • 8 ( xm;  xn ) *  ^ 1 ACn  -  ra/2  (100) 

n=H+L 


Eqs  (100)  through  (108)  represent  2N+M  equations  and  2H+M 
unknowns.  The  unknowns  are  04,  ag,  ...  a^,  Oj,  02*  •••  @n>  ti» 

r 2_%  .  .  .  and  tjj.  Each  of  the  unknowns  appear  in  each  of  the 

equations,  though  several  are  multiplied  by  zero.  Remember, 
each  of  Eqs  (100)  through  (108)  are  valid  for  a  given  set  of 
xm’s.  The  unknowns  may  be  solved  for  using  matrix  techniques. 
This  is  easier  to  see  if  the  equations  are  set  in  the  following 
form: 
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1 

1 

11 

12 

1 

1 

21 

22 

1 (2N+M)  (1  ) 


•  •it 


1 

(1  )  ( 2N+H) 

a 

1 

• 

• 

• 

a 

M 

• 

e. 

• 

• 

• 

"h 

t 

r 
.  1 

1 

(2N+M)  (2N+H)J 

T 

L  H  J 

0 

0 

8  (M+  1  ) 

g 

(2M) 

0 

0 


(109) 


Where  lmn  are  determined  from  Eqs  (100)  through  (108)  and 
represent  the  inner  products  in  Eqs  (88)  through  (91).  The 
vector  in  Eq  (109)  with  the  coefficients  a,  3,  and  r  naturally 
represents  the  unknowns.  The  coefficients  gm  represent  the 
incident  magnetic  field  in  Eq  (108)  obtained  from  Eq  (74). 
Representing  Eq  (109)  as 


[  lmn  1  (Qn^  =  I®mJ 

the  unknowns  by  can  be  found  by  the  inverse  operation  (or  its 
equivalent ), 


[<*n  J  =  UmnJ“lf«m]  (HO) 

Once  the  unknowns,  a’s,  0’s,  and  t’s,  are  determined 
numerically,  they  may  be  used  in  Eq  (61)  to  approximate  the 
currents  and  charges  as  functions  of  position.  These 
approximations  may  then  be  used  to  find  the  scattered  and 
total  electric  and  magnetic  fields  using  Eqs  (58)  and  (59), 


58 


respectively,  for  all  points  in  the  half  space  above  the  y  =  0 
plane  and  Eqs  (60)  and  (61),  res  pect  i  ve  1  y,  for  all  points 
inside  the  cavity.  The  same  approximations  used  in  the 
previous  integrals  may  be  used  to  approximate  the  integrals  in 
Eqs  (58)  and  (59).  If  the  radar  cross  section  of  the  cavity  is 
desired  instead  of  the  scattered  fields,  an  equation  needs  to 
be  developed  to  produce  it.  The  required  equation  is  derived 
in  the  next  chapter. 
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VT 


Ha d»r  Cross  Section 


In  three  dimensional  space,  the  RCS  (radar  cross  section) 
of  an  obstacle,  or  scatterer,  is  defined  as  the  area  for  which 
an  "incident  wave  contains  sufficient  power  to  produce,  by 
omnidirectional  radiation,  the  same  back-scattered  power 
density"  (4:116).  Simply  put,  RCS  is  the  cross  sectional  area 
that  would  normally  be  required  to  isotropically  scatter  the 
same  power  as  the  target  radiates  toward  the  receiver.  The 
mathematical  form  for  RCS  can  be  in  terms  of  the  incident  and 
scattered  power,  Pi  and  Ps  respectively,  as  in  (4:116): 

I JPS  I 

o  =  1  im  4irr2  -  (111) 

r-*co  j  Pi  | 

RCS  may  also  be  determined  from  the  field  quantities,  E  and  H, 
as  in  (9:  157  ): 

IF*  I  2 

<x  =  1  im  4wr2  -  (112) 

r— *00  |  F 1 1  2 

where  F  refers  to  either  E  or  H. 

Though  other  definitions  are  possible  for  RCS,  Eqs  (111) 
and  (112)  are  the  most  widely  used.  The  assumptions  that  are 
made  in  deriving  Eqs  (111)  and  (112)  are  that  the  scatterer  is 
a  three  dimensional  object  and  the  distance,  r,  from  the  object 
approaches  infinity  to  remove  the  dependence  on  range.  The 
scattered  fields  decay  as  1/r  in  the  far  field.  The  presence 
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of  r2  in  the  numerator  of  Eqs  (ill)  and  (112)  effectively 
eliminates  the  effect  of  range  on  RCS. 

For  two-dimensional  RCS  problems,  such  as  in  the  channel 
problem,  the  RCS  is  the  area  (or  width)  that  would  normally  be 
required  to  scatter  " i s o trop i c a  1 1 y" ,  in  two-dimensions,  the 
same  power  radiated  toward  the  receiver.  A  three  dimensional 
isotropic  radiator  radiates  into  a  sphere;  a  two-dimensional 
isotropic  radiator  radiates  into  a  circle.  The  fields  decay  as 
r-1/2  in  two-dimensional  problems.  Though  not  actually  a 
cross  sectional  area,  the  term  RCS  is  still  used  in  two- 
dimensions.  Perhaps  a  more  correct  term  would  be  a  scattering 
width  (9).  The  equation  for  RCS  using  the  incident  and 
scattered  fields  in  two-dimensions  becomes: 


a 


1 im  2nr 
r  a) 


|£*l 2 

ir‘|2 


(113  ) 


The  distance  from  the  scatterer  is  again  large  enough  to 
remove  the  dependence  on  range.  Notice  that  if  a  two- 
dimensional  scatterer  could  be  seen  in  three  dimensions,  it 
would  be  infinitely  long.  If  a  two-dimensional  scatterer  of 
bounded  cross  sectional  area  could  be  viewed  in  three 
dimensions  from  a  place  a  large  distance  r  away,  the 
scatterer  would  look  like  an  infinite  line. 

To  use  Eq  (113)  on  the  channel  problem,  it  must  be  noted 
that  only  the  RCS  of  the  cavity  can  be  found.  Since  the  cavity 
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is  embedded  in  an  infinite  PEC  plane,  it  would  be  impossible  to 
reach  a  distance,  r,  such  that  the  plane  would  appear  as  a 
point.  Therefore,  only  the  cavity’s  RCS  can  be  determined. 

This  is  not  at  all  bad.  Determining  the  RCS  of  an  infinite  PEC 
plane  is  nearly  useless. 

The  magnetic  field  scattered  from  the  channel  can  be  found 
using  Eq  (76)  once  the  magnetic  current  is  known.  It  appears 
as 


iwe 

Hs  (x)  =  - 

4ir 


Perhaps  more 


J  M(x  *  )G  ( x;  x  1  )  dx  * 

A 

precisely,  using  Eqs 


(70  )  through  (73  ), 


(114) 


'  1/2 
ue 

Hs (x )  =  -  -  --  M(x* )H<4 ){k[(x-x*  )2  +  y2)1/2)  dx*  (115) 

-  -  2  O 

J-i/2 

Eq  (113)  refers  to  the  limit  as  r  approaches  infinity.  Since 
Eq  (115)  is  in  (x, y)  coordinates,  it  needs  to  be  converted  to 
(r, 0)  coordinates.  If  P  is  the  angle  from  the. positive  y  axis 
measured  positively  clockwise  and  r  is  the  distance  from  the 
origin,  then  the  Hankel  function  of  the  first  kind,  order  zero 
in  Eq  (115)  becomes 

>{k[ (x-x*  )2  +  y2) l/2i 

=  H^1 ) ( k (r2s in20  -2x’rsin0  +  x’2  +  y2cos20)1/2) 


62 


o  r 

H^1  ){K[(X-X»  )2  +  y2]  1/2 1  =  H^hktr2  -  2x  •  rs  in9  +  x’2]1/2! 

(116) 

s  ince 

x  =  r  s in (9  ) 
y  =  r  cos (9 ) 

For  large  arguments,  such  as  when  r  goes  to  infinity  in  Eq 
(116),  the  Hankel  function  of  the  first  kind  order  aero  behaves 
as 


r  2  i1/2 

H(1)(2)/n^  -  exp  [  i  ( z-n/4 )  ] ,  z-*®  (117) 

°  L  W2  J 


Using  a  binomial  expansion  on  the  argument  of  Eq  (116)  and 
discarding  terms  on  the  order  of  r-1,  a  simpler  form  of  Eq 
(115)  is  possible.  In  antenna  theory  it  is  common  practice 
when  going  into  the  far  field  to  expand  the  argument  and  keep 
the  first  term  for  magnitude  purposes  and  the  first  two  terms 
for  phase  purposes.  This  is  what  has  been  done  here.  The  new 
form  of  Eq  (115)  is 


Hs(r,  9) 


iu>e 


'1/2 

r  2  -,i/2 

M(x')  -  exp  [ i  (kr-x  * ks  in9 

L  irkr  -I 

J  ~  1/2 


w/4)]  dx' 

(110) 


The  integration  in  Eq  (119)  is  in  terms  of  x’  so  several 
of  the  constants  may  be  removed  to  outside  of  the  integral 
sign.  Also,  using  Eq  (81)  to  approximate  iweH(x’)  and 
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breaking  the  integral  from  -1/2  to  1/2  into  M  integrals  each 
over  a  width  ACn,  Eq  (114)  may  be  further  simplified  as 


H*  (r,  9) 


xn+AC/2 


i  r  2  -.1/2  H 

—  -  exp[  i  (Rr-ir/4 ))  E  a 

4  L  wkr  J  n=l  n 


-  exp [- ikx ' s in ( 9 )]  dx 
xn-AC/2 

(119) 


. 


After  performing  the  integration,  the  magnitude  of  the 
scattered  field  becomes 


H  s in [k (AC/2 )s in ( 9  )  ] 

Ea  exp[-ikx  sin(9)] - 

n=ln  n  k  s in ( 9 ) 


(120) 


and  k  times  the  RCS  becomes 


ka 


4 


M  s in [k (AC/ 2  )s in (9  )] 

E  a  exp[-ikx  sin(9))  - 

n= 1  n  n  k  Sjn(9) 


2 


(121  ) 


Notice  that  when  sin(9)  =  0,  the  ratio  in  Eq  (121)  is  equal 
to  AC/2.  Once  the  coefficients  of  Eq  (81)  have  been 
determined  (i.e.  a’s,  0's,  and  t’s)  the  RCS  of  the  cavity  may 
be  generated. 
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VII.  The  Computer  Program 


A  FORTRAN  computer  program  has  been  written  to  implement 
Eqs  (100)  through  (108)  and  solve  for  the  unknown  coefficients 
in  Eq  (81).  Using  these  coefficients,  the  program  then 
evaluates  Eq  (121)  for  the  RCS  of  the  cavity.  The  fields 
scattered  from  the  cavity  or  from  both  the  PEC  plane  and  the 
cavity  are  easily  obtainable  using  Eqs  (59)  through  (62),  but 
have  not  been  plotted.  A  copy  of  the  program  is  contained  in 
Appendix  A. 

The  inputs  to  the  program  are  the  channel  depth, 
normalized  by  the  width,  the  incident  angle  of  the  incoming 
field,  and  the  operating  frequency  in  the  form  of  the  wave 
number,  k.  Here,  the  wave  number  is  also  normalized  by  the 
channel  width,  w.  The  output  of  the  program  is  the  RCS  of  the 
cavity  excluding  the  PEC  plane  at  all  reflected  angles  from 
-90°  to  +90°,  measured  from  the  +y  axis. 

Appendix  B  provides  some  of  the  output  from  the  program  in 
a  plotted  form.  The  input  angles  for  the  plane  wave  are  0^ 
equal  to  0°,  22.5°,  45°  and  67.5°.  The  wave  numbers  used  are  k 
equal  to  0.1,  1.0,  and  10.0  which  correspond  to  long 
wavelength,  intermediate  wavelength,  and  short  wavelength 
respectively.  The  normalized  depth  of  the  channel,  d,  is 
either  0.25,  1.0,  4.0,  or  8.0. 

Using  these  numbers  and  subdividing  the  channel  into 
pieces  small  enough  so  that  the  smallest  length  in  the  problem 
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(depth,  wavelength,  or  width)  is  resolved  by  at  least  12 
patches,  will  yield  systems  of  equations  (2N+M  by  2N+M)  ranging 
in  size  from  a  minimum  of  108  by  108  to  a  maximum  of  440  by 
440.  The  maximum  number  of  equations  and  unknowns  the  program 
can  reasonably  solve  is  450  due  to  storage  constraints.  A 
Gaussian  elimination  subroutine  utilizing  scaled  partial 
pivoting  and  iterative  improvement  was  used  to  solve  the  linear 
system  of  equations.  Matrix  inversion  was  tested  as  a  possible 
method  of  solution  and  found  to  be  as  accurate  as  the  Gaussian 
elimination  subroutine  (to  eight  significant  digits)  and  much 
faster  for  more  than  one  input  angle. 

Using  the  Gaussian  elimination  subroutine  was  slightly 
faster  than  the  matrix  inversion  for  one  input  angle  and  the 
output  of  both  solutions  was  comparable.  However,  a  single 
matrix  inversion  can  be  used  to  determine  the  RCS  for  many 
incident  angles.  With  more  than  one  input  plane  wave  (i.e. 

9 1  =0°,  22.5°,  45°,  and  67.5°),  it  was  more  efficient  to  use  the 

matrix  inversion  technique.  The  output  plotted  in  Appendix  B 
was  generated  using  the  matrix  inversion  technique. 

The  program  given  in  Appendix  A  can  be  .easily  modified  to 
allow  the  variables  K  and  d  to  vary  while  keeping  the  incident 
and  reflected  angles  constant.  Doing  so  requires  a  matrix 
inversion  for  each  different  k  and  d.  This  slows  the  program 
down  since  several  inversions  are  necessary  to  get  a  good  plot 
for  the  data.  The  plots  obtained  by  varying  k  and  d  for 
specific  input  and  output  angles  are  presently  more  interesting 
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than  those  plots  contained  in  Appendix  B.  Examples  of  varying 
k  and  d  are  given  in  the  next  chapter  with  a  discussion  of  the 
results  of  the  thesis. 

Table  III  gives  a  quick  list  of  the  plots  obtained  by 
running  the  program.  The  plots  generated  from  the  program's 
data  were  obtained  from  the  program  inputs  given  in  the  table. 


Table  III.  Summary  of  Inputs  to  the  Plotted  Data 


Normal ized 
Wave  Number 
k 

Normal ized 
Channel  Depth 
d 

Plane  Wave 
Incidence  Angle 

0 

i 

Matrix 

Size 

2N+M 

0. 1 

0.25 

0, 

22.5,  45, 

67.5 

288 

0.  i 

1.0 

tt 

ft 

ft 

108 

0.  i 

4.0 

It 

It 

tt 

252 

0.  i 

8.0 

ft 

tt 

N 

444 

i  .0 

0.25 

ft 

ft 

tt 

288 

1 .° 

1.0 

tt 

tt 

tt 

108 

1 . 0 

4.0 

N 

ft 

tt 

252 

1.0 

8.0 

ft 

tt 

tt 

444 

10.0 

0.25 

ft 

tt 

tt 

288 

10.0 

1.0 

tt 

tt 

tt 

144 

10.0 

4.0 

tt 

It 

N 

336 

0.2 - 10.0 

1.0 

0,  45 

- 

0.2 - 10.0 

4.  0 

0,  45 

- 

1.0 

0.25 - 8.0 

0,  ,  45 

- 

8.0 

0.25 -  8.0 

0,  45 

Total  Number 

Of  Plots  =  60 
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The  theory  and  equations  derived  earlier  allow  for 
determining  the  magnetic  or  electric  field  at  any  point  x  or 
the  "CS  at  any  bistatic  angle.  Notice  that  x  may  be  inside  the 
cavity,  on  the  cavity,  in  the  upper  half  space,  or  at  any 
desired  point  in  space.  The  plots  given,  however,  are  only  for 
the  RCS  of  the  given  input  values;  no  plots  were  generated  for 
the  fields  inside  the  cavity,  though  it  is  trivial  to  find  them 
if  desired.  Discussion  of  the  most  prominent  results  is  given 
in  the  next  chapter. 
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VIII.  Resul ts 


The  most  important  contribution  of  this  thesis  is  the 
development  of  the  coupled  set  of  Fredholm  integral  equations 
of  the  second  kind.  A  set  of  equations  was  developed  for  the 
region  above  the  cavity  and  a  set  of  equations  was  developed 
for  use  within  the  cavity.  These  equations  are  coupled 
together  by  the  aperture  of  the  cavity.  The  coupled  set  (Eqs 
(64) through  (67))  of  Fredholm  integral  equations  of  the  second 
kind  can  be  used  for  any  cavity  in  a  PEC  plane;  they  were  later 
specialized  for  use  with  a  two-dimensional  cavity  when  a  TE 
plane  wave  is  the  incident  field.  In  this  particular  case,  the 
equations  can  easily  be  solved  numerically;  or  at  least 
approximated  numerically. 

The  resulting  set  c!  equations  developed  for  the  two- 
dimensional  problem  (Eqs  (76)  through  (80))  reduced  an 
integration  over  an  infinite  plane  to  an  integration  over  an 
aperture  and  around  the  cavity  within  the  plane.  These 
specific  equations  can  be  used  with  any  shape  two-dimensional 
cavity.  The  scattered  fields  can  be  determined  by  using  Eq 
(76)  and  the  currents  and  charges  found  using  Eqs  (77)  through 
(80). 

A  specific,  rectangular,  cavity  was  introduced  and  the  RCS 
determined  from  that  shaped  cavity  for  several  wavenumbers,  k, 
and  cavity  depths,  d.  A  computer  program  was  written  and  used 
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to  solve  those  equations.  Most  of  the  data  generated  by  the 
computer  program  is  given  in  Appendix  B.  Some  of  the  most 
important  observations  of  that  data  will  be  presented  here. 

Figure  5  shows  three  plots  of  RCS  data  for  a  wavenumber  of 
10,  an  incident  angle  of  0°,  and  channel  depths  of  0.25,  1.0, 

and  4.0.  Notice  that  the  shape  of  the  curves  remains 
relatively  constant  whilv  the  peak  magnitudes  change.  With 
incident  angles  of  0°,  all  of  the  other  sets  of  k  and  d  produce 
plots  similar  to  Figure  5.  This  indicates  that  a  0°  plane  wave 
(with  any  wavelength)  incident  to  any  rectangular  cavity 
produces  a  bistatic  RCS  of  roughly  the  same  shape  given  in 
Figure  5;  only  the  magnitudes  change  with  changing  wavenumber 
and  cavity  size. 

Figure  6  shows  the  same  wavenumber  and  depths  when  the 
incident  angle  is  45°.  At  this  incident  angle,  the  shape  of 
the  curves  do  not  remain  constant.  This  inchoates  that  depth, 
wavenumber,  and  incident  angle  all  contribute  to  RCS  separately 
and  are  not  tied  specifically  together. 

Figures  7  and  8  give  plots  of  k=1.0  and  d=i.O  with 
incident  angles  of  22.5°,  45°,  and  67.5°  and  k=1.0  and  d=0.O 
with  incident  angles  of  22. 5°,  45°,  and  67.5°,  respectively. 

With  varying  incident  angles  and  a  fixed  k  and  d,  the  shape  of 
the  curves  tend  to  change  slowly.  The  relative  magnitude  of 
the  "lots  also  changes  slowly.  This  type  of  plot  occurs  for 
eery  channel  depth  when  k=0.1  and  k-i.C  When  k=i0.0  (see 
Appendix  B  for  the  plots)  the  wavelength  is  much  smaller  and  a 
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small  change  in  input  angle  changes  the  ouput  RCS  by  a  much 

I  anger  degree . 

Allowing  k  to  vary  while  holding  the  depth  and  incident 
and  scattering  angles  constant  produces  some  interesting  plots. 
Figures  9,  10,  and  11  are  examples  of  this  type  of  plot. 

Figures  9  and  10  are  monostatic  at  0°  with  d-1.0,  while  Figure 

II  is  bistatic  at  0°  and  45°.  Notice  that  Figure  11  contains 
two  curves;  one  for  an  input  angle  of  0°  (output  of  45°)  and 
one  for  an  input  angle  of  45°  (output  of  0°).  The  concept  of 
reciprocity  (21)  is  represented  quite  well  from  the  nearly 
exact  coincidence  of  the  two  plots. 

The  plot  in  Figure  9  approaches  zero  at  three  points:  all 

integer  multiples  of  n.  When  K  =  mr,  the  wavelength,  X,  is 
2/n.  With  a  depth  of  1.0  and  an  incident  angle  of  0°,  the 
monostatic  RCS  from  this  cavity  approaches  zero  because  an 
integer  number  of  half  wavelengths  exactly  fill  the  cavity. 

The  monostatic  RCS  at  45°  (Figure  10)  has  what  seems  to  be 
discontinuities  around  k  =  3.4,  k  =  6.3,  and  k=9.7.  Though  not 
shown,  these  RCS  "spikes"  occur  at  incident  angles  of  15°,  30°, 
and  60°  at  near  the  same  values  of  k,  and  expected  for  at  least 
all  angles  in  between  15°  and  60°. 

This  phenomena  was  investigated  further  by:  1)  expanding 

the  curve  around  the  k=3.4  value  (Figures  12  and  13)  and  2) 
allowing  the  depth  to  vary  with  k=3.3,  k;3.4,  and  k=3.5  (Figure 

11).  From  Figures  12  and  13,  it  appears  as  if  the  "spike"  is 
at  least  a  smooth  curve  and  may  belong  in  the  RCS  plots  (as 


opposed  to  being  a  numerical  instability).  Figure  14  shows 
that  slight  changes  in  k  and  d  ground  these  :p:'.'.es  give  large 
changes  in  the  monostatic  RCS  at  45°.  Again,  this  indicates 
that  k,  d,  and  0i  are  all  separate  variables  and  should  be 
treated  as  such. 

Holding  k  constant  and  allowing  d  to  vary  produces 
oscillating  plots  from  a  maximum  to  a  minimum  and  back  again. 
This  is  because  as  the  depth  of  the  channel  is  increased,  more 
wavelengths  (when  the  wavelength  is  held  constant)  can  fit  into 
the  channel.  An  example  of  this  type  of  plot  is  given  in 
Figure  15.  Figures  similar  to  Figures  5  through  15  are  given 


in  Append ix  B . 


k=  1 0.0  d=0.25,  1.0,  &  4.0  Inc  Angle=45' 


b 


Figure  6.  Bistatic  RCS  with  K  =  i0.0,  (1  =  0.35,  i.0,  and  4.0  with 

an  Incident  Angle  of  45°. 
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k=  1 .0  d=8.0 

Gj  =  22.5°,  45°,  &  67.5° 


Figure  8.  Bistatic  RCS  with  K=i.O  and  d^fl.O.  The  Incident 
Angles  are  22.5°,  45°.  and  67.5°. 
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Monostatic  Scatter  at  45°  with  d=1.0 


Figure  10.  Monostatic  RCS  at  45°  with  d= 
wavenumber,  k. 
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Bistatic  RCS  at  0°/45°  with  d=i,0  and  a  varying 
k.  Reciprocity  can  easily  been  seen  from  these  two 
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d=1.0  Monostatic  Scatter  at  45 


2.7  2.9  3.1  33  35  37  3® 

k  (amal  Incramants) 


Figure  12.  Monostatic  RCS  at  45°  with  d=1.0  and  a  varying 
wavenumber,  k;  Small  increments  of  k  were  taken  around  the 
spike  in  Figure  10. 
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d=1.5  Monostatic  Scatter  at  45° 


Figure  13.  Monostatic 
wavenumber,  k;  Small 
spike  in  Figure  10. 


RCS  at  45°  with  d  =  1.5  and 
increments  of  k  were  taken 


a  varying 
around  the 
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b 

X 


Monostotic  Scatter  at  45° 


Varying  d  and  l«a  aping  k  oonakmt 


Figure  14.  Monostatic  RCS  at  45°  with  k  =  3.3,  3.4,  and  3.5  and 

a  varying  depth,  d;  Small  increments  of  d  were  taken  for  the 
three  k  values  around  the  spike  in  Figure  10 . 
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Monostotic  Scatter  at  O^With  k=1.0 


Choral  Dapth,  d 


Figure  15.  Monostatic  RCS  at  0°  with  a  fixed  K=i.O  and  a 
varying  depth,  d. 
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IX ■  Cone  1  us  iona  and  Re commend at  ion a 

The  most  important  contribution  of  the  thesis  is  the 
development  of  the  coupled  set  of  Fredholm  integral  equations 
of  the  second  Rind  for  the  two-dimensional  cavity.  The 
integral  equations  derived  for  finding  both  the  electric  and 
magnetic  fields  and  the  surface  currents  and  charges  are  exact. 
This  thesis  shows  that  using  the  half  space  Green’s  function  to 
derive  the  equations  above  the  plane  in  Figure  1,  the  free 
space  Green’s  function  to  derive  the  equations  used  inside  the 
cavity,  and  then  coupling  these  equations  through  the  aperture 
is  an  accurate  way  to  determine  the  fields  scattered  from  the 
channel.  It  is  concluded  that  the  RCSs  plotted  in  Figures  5 
through  15  and  Appendix  B  are  reasonably  accurate.  They  are 
not  exact  in  that  some  approximations  were  necessary  to  solve 
the  equations  numerically. 

The  method  used  to  calculate  the  fields,  currents,  and 
charges  is  an  approximation.  As  the  discretization  parameter, 
AC,  approaches  zero,  the  approximation  should  better  represent 
the  exact  solution.  However,  with  more  segments,  the  number  of 
unknowns  and  equations  also  increases  and  therefore  so  does  the 
computer  time  required  to  solve  the  system  of  equations.  As 
the  number  of  equations  increases,  so  too  does  round-off  error 
in  the  computer.  The  accuracy  that  can  be  obtained  using 
moment  methods  is  therefore  limited  by  round-off  error. 

In  this  thesis,  the  smallest  length  in  the  problem 
(wavelength,  channel  depth,  or  width)  was  broken  up  into  at 


94 


least  12  segments.  This  seems  to  give  reasonable  accuracy  at  a 
reasonable  cost  in  computer  time.  Test  cases  were  applied  for 
several  different  segments  per  smallest  length.  (Actually,  4, 

5,  6,  8,  10,  12,  14,  16,  20,  24,  and  32  segments  per  smallest 

length  were  used. )  All  of  the  cases  above  8  segments  per 
smallest  length  gave  RCS  plots  of  approximately  the  same 
magnitude  and  shape  for  all  bistatic  angles.  The  time  required 
when  using  12  segments  per  smallest  length  was  only  a  -*ew 
minutes  longer  than  when  using  10  segments.  Though,  even  at  32 
segments  per  smallest  length,  the  RCS  plots  seemed  to  be  still 
converging  on  an  "exact"  plot  (remember,  numerical  solutions 
are  limited  by  round-off  error.  )  The  theory  and  methods  used 
here  allow  for  more  accuracy,  limited  by  round-off  error,  if 
desired.  The  user  must  be  willing  to  give  up  computer  time  in 
order  to  achieve  this  accuracy. 

Appendix  C  gives  a  description  of  a  target  designed  to 
obtain  RCS  measurements  of  the  rectangular  cavity.  Using  this 
target  and  the  method  described  in  Appendix  C,  laboratory 
measurements  can  be  obtained  to  check  the  predicted  RCS  values 
given  in  this  thesis. 

There  are  several  possibilities  for  further  analysis  for 
this  problem.  These  include: 

1.  Making  the  measurements  on  the  target  described  in  Appendix 
C  to  determine  the  accuracy  of  the  method  described  therein. 

2.  A  selection  of  different  shaped  cavities  than  the 
rectangular  one  used  here.  This  would  require  reworking  the 
approximate  solutions  and  new  computer  programs. 
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3.  A  TM  (transverse  magnetic)  incident  field  instead  of  the  TE 
case  used  here.  New  equations  may  have  to  be  derived  using 
Dinchlet  boundary  conditions  and  an  equation  similar  to  Eq 
(58)  could  be  used  in  place  of  Eq  (60)  and  an  equation  similar 
to  Eq  (61)  used  instead  of  Eq  (59)  to  get  equations  similar  to 
Eqs  (64)  and  (65)  respectively. 

4.  A  cavity  filled  with  a  dielectric.  This  changes  the 
problem  drastically.  The  free  space  green’s  function  inside  of 
the  cavity  will  have  to  be  evaluated  to  determine  the  affects 
of  filling  the  cavity  with  a  dielectric. 

5.  Change  the  channel  to  be  an  open  channel,  meaning  that  the 
depth  is  infinite.  New  equations  would  need  to  be  derived  for 
the  inside  of  th«  cavity  in  this  case.  Perhaps  a  new  Green’s 
function  can  be  explored  to  eliminate  integrating  along  the  now 
infinitely  deep  channel  sides. 

6.  Let  the  cavity  be  three  dimensional.  That  is,  let  the 
cavity  have  a  finite  length  as  well  as  the  finite  width,  w,  and 
depth,  d.  This  gets  closer  to  an  open  ended  flanged  waveguide. 


Another  contribution  from  this  thesis  is  that  it  may  now 
be  possible  to  easily  determine  the  RCS  from  infinite 
cylinders,  or  other  two-dimensional  objects,  containing  abrupt 
"cavity-1  ike"  variations  in  an  otherwise  smooth  surface.  If 
such  an  object  had  a  trench  along  the  side,  such  as  the 
rectangular  cavity  explored  in  this  thesis,  or  some  other 
cavity,  it  should  be  possible  to  find  the  RCS  of  such  a 
cylinder  by  adding  the  fields  scattered  from  a  perfect  cylinder 
to  those  scattered  from  the  perturbing  cavity  and  using  the 
result  as  the  scattered  field.  Such  an  approximation,  based  on 
the  cavity  solution  obtained  herein,  would  better  represent  the 
true  solution  when  the  curvature  near  the  cavity  is  small. 
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Figure  16  shows  the  proposed  method  to  approximate  the 


scattered  RCS  from  a  cylinder  with  a  small  rectangular  cavity. 
It  should  be  investigated  to  ascertain  the  feasibility  of  the 
method . 


(a)  Fields  Scattered  from  a  Whole  Cylinder  Without  the  Cavity 
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(b)  Find  the  Fields  Scattered  from  the  Cavity 


(c)  Add  the  Fields  from  (a)  and  (b)  to  Calculate  the  RCS 


Figure  16.  Method  for  Finding  the  RCS  cf  an  Infinite  Cylinder 
With  a  Rectangular  Cavity  by  Adding  the  Fields  Scattered  From 
the  Perfect  Cylinder  (a)  to  Those  Scattered  From  the  Cavity  (b) 
to  Get  the  Total  Field  (c). 
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The  above  recommendations  for  further  study  are  important 
to  fully  evaluate  the  scattered  EM  fields  from  objects 
containing  cavities.  The  information  contained  in  this  thesis 
is  a  continuation  toward  that  goal. 

In  the  preceding  pages,  a  set  of  Fredholm  integral 
equations  of  the  second  kind  using  Neumann  boundary  conditions 
(the  half  space  Green’s  function,  G(x;x*))  were  derived  for  the 
half  space  above  the  cavity.  These  equations  were  then  coupled 
to  a  corresponding  set  of  integral  equations  for  inside  the 
cavity  and  the  coupled  set  was  used  to  find  surface  currents 
and  charges  that  could  be  used  to  find  the  scattered  fields  or 
the  RCS  of  the  cavity.  A  specific,  rectangular,  cavity  was 
introduced  as  a  test  case  and  the  RCS  calculated  using  a 
Fortran  computer  program  specifically  written  for  this  task. 

The  output  of  that  program  was  addressed  and  presented. 

Finally,  specific  areas  for  future  analysis  were  given.  Though 
the  analysis,  study,  and  experimentation  continues,  this  thesis 
has  completed. 
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The 

written 


Appendix  Ai  Computer  Program 

following  listing  is  the  FORTRAN  computer  program 
to  solve  the  test  problem  within  the  thesis. 
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C  This  is  the  program  for  finding  the  RCS  of  a  rectangular 

C  channel  of  normalized  depth  (dw),  unity  width,  and  for  a 

C  wave  number  (wk)  read  as  input.  This  version  runs  on  the 

C  AFIT  CSC  computer  system  with  IMSL  commands.  The  program 

C  must  be  linked  with  the  IMSL  library.  To  do  this  on  the 

c  CSC 

C 

C  THE  APPROPRIATE  LINE  STATEMENT  IS: 

C  LINK  FILENAME. OBJ,  IMSLiDIR:  IMSL/LIB 

C 

C  This  version  uses  the  matrix  inverse  routine. 

C  It  reads  the  input  file  "DINP.DAT"  for  the  number  of 

c  different  K’s  (wk )  and  d’s  (dw)  to  be  used.  This  number  is 

c  called  nin.  The  program  then  runs  nin  times  and  prints  out 
c  the  RCS  data  in  10  significant  digits.  The  number  of 
c  segments  that  the  smallest  measurement  is  to  be  broken  into 

c  is  read  from  an  input  file  as  well  as  the  incoming 

c  planewave  angles. 

c  This  version  checks  for  1/delc  as  an  integer.  As  long 

c  as  npp  is  divisible  by  4,  d  is  a  multiple  of  .25  and 

c  K  =  1 0 ,  or  less  than  2pi,  this  will  work. 

c 

C  Some  of  the  main  variables  are: 
c  J: imaginary  operator  (squareroot  of  -1) 
c  dw=channel  depth  normalized  to  the  width 
c  wK:wave  number  normalized  to  the  width 
c  delc=width  of  each  line  segment  after  discretization 
c  RCS();Radar  Cross  Section  data  that  is  output 

C  XM(l(or  2),  N)=The  midpoint  of  the  Nth  linesegment,  a  1 

c  indicates  the  x  point  and  a  2  indicates  the  y  point 

c  CLNM (m,  n  )  =  the  matrix  holding  the  data  for  the  linear 

c  operations  on  the  currents  and  charges.  (see  the  Thesis  c 

noted  above  Eq.  (109)  for  details) 

cout()=the  current  and  charge  coefficients 
c  g()  =  the  incident  fields  data  (see  Eq.  190  in  the  thesis) 
c 
c 

IMPLICIT  DOUBLE  PRECISION  A,  B,  D,  E,  F,  H,  0,  Q,  R,  T,  U,  V,  W,  X,  Y,  Z 

IMPLICIT  COMPLEX*  16  C,  J,  G,  P,  S 
DOUBLE  PRECISION  CONST 
CHARACTER* 10  FLNAME 

DIMENSION  G  ( 1 :  500  ),  CLNM ( 1 :  500,  1:  500  ),  COUT(  1 :  50  0  ) 

DIMENSION  CWKS  (  125250  ),  IWKS  (500  ),  RCS  (  1:  10,  -90:  90  ) 

COMMON  J,  XM(  1:  2,  1:  500  ),  WK,  DELC 
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c 

C  Set  some  of  the  preliminary  constants 
c 

j:  (0  .  DO,  1  .  DO  ) 

ONE= 1 . DO 

CZERO: (0. DO, 0. DO ) 

PI :3 . 14159265359 DO 
c 

c  Open  files  for  input  (2  and  4)  and  output  (5) 
c  unit  5  holds  the  real  and  imaginary  parts  of  the  current 
c  across  the  aperture  for  the  first  incident  plane  wave 
c 

OPEN  (UNIT: 2,  FILE= 'DINP. DAT',  STATUS: 'UNKNOWN'  ) 

OPEN  (UNIT  =  4,  FILE: • DIN2 . DAT*,  STATUS: 'UNKNOWN'  ) 

OPEN  (UNIT:5,  FILE: 'CRNT.DAT',  STATUS: 'UNKNOWN'  ) 
c 

c  NIN  is  the  number  of  times  to  read  the  input  file  for  sets  of 
c  K  and  d  (which  are  wk  and  dw)  and  output  files  to  be  unit  #3 
c  NPP  is  the  number  of  points  per  smallest  measurement 
c 

READ  (2,  899  )  NIN 
READ  (2,899  )  NPP 
c 

c  RFRST  is  the  first  incidence  angle,  RLST  is  the  last,  and 

c  RSTEP  is  the  step  size  to  get  from  RFRST  to  RLST 


c 

898 

899 

C 

c 

c 

c 


C 

C 

c 


C 

C 

c 

c 


c 


READ  (4,  «  )  RFRST.  RLST,  RSTEP 
FORMAT (FI  0. 3 ) 

FORMAT (14) 

IMAX  is  the  integer  to  be  used  latter  to  store  the  data 
generated  as  RCS  data 

IMAX  = I  DINT ( (RLST-RFRST )/RSTEP+ 1 .5) 

DO  650  ICNT=  1,  NIN,  1 
READ  (2,  850  )  WK 
READ  (2,851)  DW 
READ  (2,855)  FLNAME 

OPEN  (UNIT: 3,  F I LE: FLNAME,  STATUS: 'UNKNOWN'  ) 

WK  =  DBLE (WK ) 

DW: DBLE ( DW ) 

WAVE  is  the  wavelength 

RMIN  is  the  smallest  of  WAVE,  DW  (depth),  and  1  (width) 
WAVE : 2 . D0«  PI/WK 
RMIN:  DMIN1  (WAVE,  DW,  ONE) 

Find  the  discretization  size,  DELC 

and  then  make  sure  that  1/delc  is  an  integer 

DELC:RMIN/DBLE (FLOAT (NPP )  ) 

DELC : 1 . DO/ DBLE ( FLOAT ( I  DINT ( 1 . D0/DELC+ 0 . 9997 5 )  )  ) 
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C  NMAX  will  be  the  matrix  size  (nmax  by  nmax ) .  It  is  the  total 
c  number  of  line  segments  needed  to  go  around  the  cavity  twice 
c  and  across  the  aperture  one  extra  time. 

c  The  information  is  written  out  just  to  keep  track  how  the 
c  program  is  progressing, 
c 


NMAX: I  DINT ( (4. D0«DW+5. DO )/DELC  +0 . 5D0 ) 

WRITE  ( » ,  854  )  WK,  DW,  DELC,  NMAX,  FLNAME 
WRITE  (3,854)  WX,  DW,  DELC,  NMAX,  FLNAME 

854  FORMAT(  IX,  3E14.  6,  2X,  13,  AiO  ) 

850  FORMAT (F9. 6) 

851  FORMAT ( F4 .  1  ) 

853  FORMAT (13) 

855  FORMAT (AIO) 

C 

C  FILL  THE  POSITION  MATRIX  XM( ) 

C 

C  THE  "+0.5"  INSURES  THAT  1.9999  DOES  NOT  GET  TRUNCATED 

C 

c  The  MONE  through  MEIGHT  are  the  last  line  segments  (dele’s) 
c  for  each  major  segment,  across  the  aperture  for  M,  across  the 
c  aperture  of  K,  down  S-  for  X,  ... 
c 

MONE: IDINT ( 1 . D0/DELC+  0 . 5D0 ) 

MTWO : 2  *  MONE 

MTHR :MTWO+ I  DINT ( DW/DELC+  0 . 5D0 ) 

MFOUR : I  DINT ( (3 . DO+DWJ/DELC+O . 5D0 ) 

MF I VE :MFOUR+ I  DINT ( DW/DELC+  0 . 5 DO ) 

MS  I X :MF I VE+MONE 

MSEVN:MS I X+ I  DINT ( DW/DELC+  0 . 5D0 ) 

ME I GHT : MS E VN+ MONE 

The  first  section  is  across  the  aperture  (for  Mag  Crnt) 

DO  10,  I  =  1,  MONE,  1 

XM( 1,  I ) : -0 . 5D0+ ( DBLE ( FLOAT ( I ) )-0 . 5D0 )*DELC 
XM ( 2,  I  )  =  0 . 0D0 
0  CONTINUE 

The  second  section  is  across  the  aperture  for  the  k  vector 
DO  12,  I  :MONE+ 1,  MTWO,  1 

XM ( 1 ,  I ) : 0 . 5 DO -DELC*  ( DBLE ( FLOAT ( I -MONE ) )-0.5D0) 

XM (  2,  I  )  :0 . 0D0 
2  CONTINUE 

Section  3  is  for  S- 

DO  14,  I  :MTWO+ 1,  MTHR,  1 

XM ( 2,  I )  =  - DELC* (DBLE (FLOAT (I -MTWO) )-0.5D0) 

XM ( 1 ,  I ) :-0 . 5D0 
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CONTINUE 


14 

C 

C  Section  4  is  for  sd 
C 

DO  16,  IrMTHR+1,  MFOUR,  1 

XM ( 1 ,  I )  =  -0 . 5D0+DELC* ( DBLE ( FLOAT ( I -MTHR ) )-0. 5D0 ) 
XM  (  2,  I  )  =  -DW 
16  CONTINUE 

C 

C  Section  5  is  for  s+ 

C 

DO  10,  I  =MFOUR*  1,  MFIVE,  1 

XM( 2,  I )  =  -DW+  DELC*  ( DBLE ( FLOAT ( I -MFOUR ) )-0. 5D0 ) 

XM(  1,  1)^0. 5 DO 
18  CONTINUE 

C 

C  Section  6  is  for  A  for  the  sigma  equations 
C 

DO  22,  I  =MF I VE+  1 ,  MSIX,  1 

XH ( 1 ,  I  )=0. 5D0-DELC*  (DBLE (FLOAT ( I -MFIVE )  )-0.5D0) 
XM  (  2,  I  )  =  0 . 0D0 
22  CONTINUE 

C 

C  Section  7  is  for  s- 
C 

DO  24,  I  =MS I X+ 1 ,  MSEVN,  t 

XM ( 2,  I  )  =  -DELC*  ( DBLE (FLOAT ( I -MSIX) )-0.5D0) 

XM ( 1 , I )=-0.5D0 
24  CONTINUE 

C 

C  Section  0  is  for  sd 
C 

DO  2  6,  I  :MSEVN+  1,  MEIGHT,  1 

XM( l,  I  )  =  -0 . 5D0+DELC*  (DBLE (FLOAT ( I-MSEVN)  )-0 . 5 DO ) 
XM (  2,  I  )  =  -DV 
26  CONTINUE 

C 

C  Section  9  is  for  s+ 

C 

DO  20,  I=MEIGHT+1,  NMAX,  1 

XM ( 2,  I  ) = - DW+  DELC »  ( DBLE ( FLOAT ( I -ME I GHT ) J-0.5D0) 

XM ( 1 ,  I  )=0. 5 DO 
20  CONTINUE 

C 

C  The  position  matrix  is  now  filled 
C 
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C  Next  we  must  fill  the  CLNM( )  matrix. 

C 

C  Set  constants 
C 

EULER= . 577 21 566490 15D0 
CONST: 2 . DO/WK/P I 
ETA=WK*  DELC/4 . DO 

SP=-CONST*(  DLOG (ETA ) - 1 . DO  ) *  ETA 
SP  =  SP+(  J/WK  -  CONST* EULER  ) *  ETA 
ETA3  =  ETA*  ETA*  ETA 

SP : SP+ ETA3 » CONST* (  DLOG (ETA)  -  i . DO/3 . DO  )/3 . DO 
SP  =  SP-ETA3* (  CONST* (1. DO -EULER)  +  J/WK  )/3.D0 
C3:-WK*WK 
C 

C  The  M  is  down  the  side  of  the  matrix  (row)  and  refers  to  the 
C  "OBSERVATION"  point.  The  "N"  is  across  the  top  of  the  matrix 
C  (column)  and  is  the  "SOURCE"  point, 
c 

c  In  appropriate  cases,  the  program  must  check  for  self- 
C  patches,  i.e.  those  cases  when  the  midpoint  of  the 
c  observation  point  is  equal  to  the  midpoint  of  the  source 
c  point. 

C 

C  THE  FIRST  MAJOR  SECTION  HAS  M  ON  A (ALPHA) 

C 

C  EQ  A  (in  the  derivation  notes) 

C 

DO  90,  M=  1,  MONE,  1 
C 

C  N  on  A(alpha) 

C 

DO  40,  N=  1,  MONE,  1 

IF  (N  .EQ.  M)  THEN 

CLNM  (M,  N):  (0. 5D0,  0  .  DO  ) 

GOTO  40 
END  IF 

CLNM (M,  N) rCZERO 
40  CONTINUE 

C 

C  N  ON  A (BETA) 

C 

DO  50,  N  =  MONE*  1,  MTWO,  1 

IF  ( DAB 3 ( XM ( 1 ,  M ) - XM ( i  ,  N )  )  . LT .  DELC/3 . DO )  THEN 
CLNM(M,  N)  :C3*SP 
GOTO  50 
END  IF 

CLNM  (M,  N  )  :  C 3  *  PH  1 0  ( XM  (  1 ,  M  ) ,  XM  (  2,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )*DELC 
50  CONTINUE 

C 

C  N  ON  S- (BETA) 

C 

DO  60,  N:MTWO+ 1,  MTHR,  1 
CLNM (M,  N)  cCZERO 
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CONTINUE 


60 
C 

C  N  ON  SD  ( BETA )  THROUGH  S+ ( BETA ) 

C 

DO  65,  N=MTHR+  1,  MFOUR,  1 

CLNM  (H,  N )  =  -  C  3  »  PH  1 0  ( XM  ( 1 ,  M ),  XM(  2,  M ),  XM(  1 ,  N  ),  XH  ( 2,  N )  )«DELC 

65  CONTINUE 

DO  66,  N  =  MF0UR+  1,  HFIVE,  1 
CLNM(M,  N)  -CZERO 

66  CONTINUE 
C 

C  N  ON  A (GAMMA ) 

C 


DO  70,  N  =  MF  I VE+ 1 ,  MS  I X,  i 

IF  (DABS (XM( 1,  M)-XM( 1,  N) )  . LT .  DELC/3 . DO )  THEN 
CLNM (M,  N )  =  CZERO 
GOTO  70 
ENDIF 

CLNM  (M,  N  )  =  PH  1 1 X  ( XM  (  i ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )  *  DELC 
70  CONTINUE 

C 

C  N  ON  S- (GAMMA)  THROUGH  S-f  (GAMMA) 

C 

DO  80,  N-MS I X+ i  ,  MSEVN,  1 

CLNM  (M,  N ) =  PH 1 1 X ( XM  ( 1 ,  M ) ,  XM  ( 2,  M ),  XM  ( 1 ,  N  ) ,  XM  (  2,  N  )  )  *  DELC 
80  CONTINUE 

DO  85,  N  =  MSEVN+ 1 ,  NMAX,  1 

CLNM  (M,  N  )  =  PH  1 1 X  ( XM  ( 1 ,  M  ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )»DELC 
85  CONTINUE 

90  CONTINUE 
C 

C  Thus  ends  the  section  for  M  on  A ( a  1 pha ) - - the  next  major  c 
section  is  for  M  on  A(beta) 

C 

C  Eq  b  divided  by  2  to  get  the  1/2  along  the  diagonal 
C 

DO  130,  M  =  MONE+ 1,  MTWO,  1 

C 

C  N  on  A(alpha) 

C 

DO  100  N=  1,  MONE,  1 

IF  ( DABS ( XM ( 1 ,  M ) -XM ( 1 ,  N ) )  . LT .  DELC/3 . DO )  THEN 
CLNM (M,  N)  =-SP 
GOTO  100 
ENDIF 

CLNM  (M,  N)  =  -PHI0(XM(1,  M ),  XM(  2,  M ),  XM  ( 1 ,  N  ) ,  XM  ( 2,  N )  )»DELC 
100  CONTINUE 

C 

C  N  on  A(beta ) 

C 

DO  110,  N-MONE+ 1 ,  MTWO,  1 
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IF  (N  .EQ.  M)  THEN 

CLNM(K,  N)  =  (0.5 DO,  0  .  DO  ) 
GOTO  110 
END  IF 

CLNM (M,  N)  =CZERO 
110  CONTINUE 

C 

C  N  on  all  of  the  rest  of  the  segments 
C 

DO  120,  N  =  MTWO+l,  NMAX,  1 
CLNM(M,  N)=CZERO 
120  CONTINUE 

130  CONTINUE 


C 

C  This  next  major  section  takes  care  of  the  obsrvtn  point,  M, 

C  on  s-(beta)  through  s+(beta)  for  the  source  point,  N,  on 
c  A(alpha)  and  A(beta) 

C 

C 

C  The  next  major  section  is  for  M  on  s-(beta)  and  N  on  all 
c  surfaces  except  for  A(alpha)  and  A(beta) 

C 

C  Eq  c 
C 

DO  210,  M-MTWO+ 1 ,  MTHR,  1 
C 

C  N  on  A ( a  1 pha ) 

C 

DO  140,  Nrl.MONE,  1 

CLNM  (M,  N  )  =  PH  1 0  ( XM  (  1 ,  M ) ,  XM  (  2 ,  M  ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )  •  DELC 
lao  CONTINUE 

C 

C  N  on  A ( beta  ) 

C 

DO  150,  N=MONE+ 1,  MTWO,  1 

CLNM  (M,  N  )  :  PHI  1 Y  ( XM  ( 1 ,  M  ) ,  XM  (  2 ,  M ) ,  XM  (  1 ,  N  ) ,  XM  (  2,  N  )  )  *  DELC 
150  CONTINUE 

C 

C  N  on  s-(beta) 

C 

DO  170,  N-MTWO+ 1 ,  MTHR,  1 
IF  (N  .EQ.  M)  THEN 

CLNM (M,  N)  =  (0 . 5D0,  0  .  DO  ) 

GOTO  170 
ENDIF 

CLNM (M,  N ) =CZERO 
170  CONTINUE 

C 

C  N  on  sd { beta  ) 

C 
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DO  180,  N=MTHR+ 1 ,  MFOUR,  1 
CLNM  (M,  N )  =  -  PH  1 1 Y  ( XM  (  1 ,  M ) ,  XM  (  2,  H ),  XM  ( 1 , 
180  CONTINUE 

C 

C  N  on  s  + (beta ) 

C 

DO  190,  N=MFOUR+l,  MFIVE,  1 
CLNM  (H,  N )  =  PH 1 1 X  ( XM  ( 1 ,  M ) ,  XM  (  2,  M ) ,  XM  ( 1 , 
190  CONTINUE 

C 

C  N  on  A(gamma)  through  s+ (gamma) 

C 

DO  200,  N=MFIVE+1,  NMAX,  1 
CLNM(M,  N)=CZER0 
200  CONTINUE 

210  CONTINUE 


C 

C 

C 

C 

C 

C 

C 

C 


The  next  major  section  is  for  M  on  sd(beta) 
Eq  d 

DO  2  60,  M=MTHR+  1,  MFOUR,  1 
N  on  A ( a  1 pha ) 


DO  141,  N=  1,  MONE,  1 

CLNM  (M,  N )  =  -  PH  1 0  ( XM  ( 1 ,  M ) ,  XM  ( 2 ,  M ) ,  XM  ( 1 ,  N 
141  CONTINUE 

C 

C  N  on  A (BETA) 

C 


DO  151,  N  =  MONE+  1,  MTWO,  1 
CLNM(M,  N )  =  -  PHI  1 Y  ( XM ( 1 ,  M),  XM  (  2,  M ) ,  XM(  1 , 
151  CONTINUE 


C 

C  N  on  s-(beta) 

C 

DO  220,  N  =  MTWO+  1,  MTHR,  1 
CLNM  (M,  N  )  =  PHI  1 X  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N 
220  CONTINUE 

C 

C  N  on  sd ( beta  ) 

C 

DO  230,  N  =  MTHR+  1,  MFOUR,  1 
IF  (N  .EQ.  M)  THEN 

CLNM(M,  N)-  (-0. 5D0,  0  .  DO  ) 

GOTO  230 
END  IF 

CLNM(M,  N)  =CZERO 
230  CONTINUE 


),  XM(  2,  N )  )  *  DELC 


),  XM (2,  N)  ) m  DELC 


),  XM  ( 2,  N)  )  »  DELC 


),  XM  ( 2,  N)  ) «  DELC 


),  XM(  2,  N)  )  *  DELC 
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c 

C  N  on  s  + (beta  ) 

C 

DO  240,  N  =  MFOUR+ 1,  MFIVE,  1 

CLNM  (M,  H )  =  -  PH  I  i  X  ( XM  ( 1 ,  M  ) ,  XM(2,H),XM(i,N),XM(2,N)  )«DELC 
240  CONTINUE 

C 

C  N  on  A(gamma)  through  s  + (gamma) 

C 

DO  250,  N=MFIVE+ 1,  NMAX,  1 
CLNH(M,  N ) =CZERO 
250  CONTINUE 

260  CONTINUE 
C 

C  The  next  major  section  is  for  H  on  s+(beta) 

C 

C  Eq  e 

C 

DO  310,  M=MFOUR+l,  HFIVE,  1 

C 

C  N  on  A(alpha) 

C 

DO  142,  N  =  1 ,  MONE,  1 

CLNM  (M,  N  )  =  -  PH  1 0  ( XM  ( 1 ,  M ) ,  XM  (  2,  M ) ,  XM  ( 1 ,  N  ) ,  XM  ( 2 ,  N  )  )  *  DELC 
142  CONTINUE 

C 

C  N  on  A (BETA) 

C 

DO  152,  N=MONE+  1,  MTWO,  1 

CLNM (M,  N):-PHI1Y(XM(1,  M),XM(2,M),XM(1,N),XM(2,N)  )«DELC 
152  CONTINUE 

C 

C  N  on  s- (beta ) 

C 

DO  270,  N  =  MTW0+  1,  MTHR,  1 

CLNM(M,  N  )  =  PH  1 1 X  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )*DELC 
270  CONTINUE 

C 

C  N  on  sd(beta) 

C 

DO  280,  N  =  MTHR+ 1,  MFOUR,  1 

CLNM  (M,  N )  =  PH  1 1 Y  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N ) ,  XM  (  2 ,  N  )  )»DELC 
280  CONTINUE 

C 

C  N  on  s+(beta) 

C 

DO  290,  N  =  MF0UR+ 1,  MFIVE,  1 
CLNM  (M,  N  )  =  CZERO 

IF  (N  .EQ.  M)  CLNM  (M,  N  )  =  (  -  0 . 5D0,  0  .  DO  ) 

290  CONTINUE 
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c 

C  N  on  A(gamma)  through  s+(gamma) 

C 

DO  300,  N=MFIVE+1,  NMAX,  1 
CLKM(M,  N)=CZERO 
300  CONTINUE 

310  CONTINUE 

C 

C  All  of  the  elements  for  the  observation  point  on  the  "beta" 
c  segments  are  now  filled. 

C 

C  The  next  major  section  is  to  have  M  on  all  the  "gamma" 
c  segments  and  the  source  point,  N,  on  the  beta  segments. 

C 

DO  330,  M=MFIVE+ 1,  NMAX,  1 
DO  320  N=MONE+  1,  MFIVE,  i 
CLNM(M,  N) =CZERO 
320  CONTINUE 

330  CONTINUE 

C 

C  Now  we  again  systematically  fill  M  on  A(gamma) 

C 

C  Eq  f 
C 

DO  390,  M=MF  I VE+ 1 ,  MS  I X,  i 
C 

C  N  on  A (ALPHA) 

C 

DO  340,  N=l,  MONE,  1 

I F ( DABS ( XM ( 1 ,  M ) - XM ( 1 ,  N ) )  . LT .  DELC/3 . DO )  THEN 
CLNM (M,  N ) =CZERO 
GOTO  340 
END  IF 

CLNM(M,  N  )  =  PH  1 1  X  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )*DELC 
340  CONTINUE 

C 

C  N  on  s-(beta)  and  s+(beta) 

C 

DO  342,  N  =  MTW0+  ! ,  MTHR,  i 

CLNM  (M,  N)  =  -C3*PHI0  (XM(  1,  M),  XM(  2,  M),  XM(  1,  N ),  XM  ( 2,  N  )  )»DELC 
342  CONTINUE 

DO  344,  N  =  MFOUR+l,  MFIVE,  1 

CLNM(M,  N  )  =  C  3  *  PH  1 0  ( XM  ( 1 ,  M ) ,  XM(2,M),XM(i,N),XM(2,N)  )»DELC 
344  CONTINUE 

C 

C  N  on  A(GAMMA) 

C 

DO  350,  N  =  MFI VE+ 1 ,  MS  I X,  1 
CLNM(M,  N )  =CZERO 

IF  (N  .EQ.  M)  CLNM(M,  N)  =  (  -  0 . 5D0,  0  .  DO  ) 

350  CONTINUE 
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c 

C  N  on  s- (gamma) 

C 

DO  360,  N=MS  I X+ 1 ,  MSEVN,  i 

CLNM(M,  N )  =  -  PH  1 1 Y  ( XM  ( 1 ,  M ) ,  XM  ( 2 ,  M ) ,  XM  (  i ,  N  ) ,  XM  ( 2 ,  N )  )»DELC 
360  CONTINUE 

C 

C  N  on  sd(gamma) 

C 

DO  370,  N  =  MSEVN+ 1,  MEIGHT,  1 

CLNM  (M,  N  )  =  -  PH  1 1 Y  ( XM  ( 1 ,  M  ) ,  XM(2,M),XM(i,N),XM(2,N)  )«DELC 
370  CONTINUE 

C 

C  N  on  s  +  ( gamma ) 

C 

DO  300,  N=MEIGHT+1,  NMAX,  i 

CLNM(M,  N )  =  -  PH  1 1 Y  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  ( 2 ,  N  )  )«DELC 
380  CONTINUE 

390  CONTINUE 
C 

C  The  next  section  is  for  M  on  s-(gamma) 

C 

C  Eq  g 
C 

DO  450,  M  =  MSIX+1,  MSEVN,  i 
C 

C  N  on  A (ALPHA) 

C 

DO  400,  N=  1,  MONE,  i 

CLNM  (M,  N  )  =  PH  1 1 Y  ( XM  ( 1 ,  M ) ,  XM  ( 2 ,  M ) ,  XM  (  i  ,  N  ) ,  XM  (  2 ,  N  )  )*DELC 
400  CONTINUE 

C 

C  N  on  A (BETA) 

C 

DO  402,  N=MONE+ 1 ,  MTWO,  1 

CLNM  (M,  N)=C3#PHI0(  XM  ( 1 ,  M ) ,  XM(2,  M),  XM(  1,  N  ),  XM(  2,  N  )  )  *  DELC 
402  CONTINUE 

C 

C  N  on  sd(beta) 

C 

DO  404,  N=MTHR+ 1 ,  MFOUR,  1 

CLNM(M,  N)  =  -C3*PHI0(  XM  ( 1 ,  M ) ,  XM  ( 2 ,  M ) ,  XM  ( 1 ,  N ) ,  XM  (  2 ,  N  )  )*DELC 
404  CONTINUE 

C 

C  N  on  A ( GAMMA ) 

C 

DO  410,  N  =  MFI  VE+ i ,  MS  I  X,  1 

CLNM  (M,  N  )  =  PH  1 1 X  ( XM  ( 1 ,  M  ) ,  XM  ( 2,  M ) ,  XM  ( 1 ,  N  ) ,  XM  ( 2,  N )  )»DELC 
410  CONTINUE 
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c 

C  N  on  s-(gamraa) 

C 

DO  420,  N=MSIX+ 1,  MSEVN,  i 
CLNM(M,  N)=CZERO 

IF  (N  .EQ.  M)  CLNM(H,  H)  =  CO.  5D0,  0.  DO  ) 

420  CONTINUE 

C 

C  N  on  sd(gamma) 

C 

DO  430,  N  =  MSEVN+ 1,  MEIGHT,  1 

CLNM(M,  N  )  =  PH  1 1 X  ( XM  ( 1 ,  M ) ,  XM  (  2,  H ) ,  XM  ( 1 ,  N ) ,  XM(  2,  N  )  )*DELC 

430  CONTINUE 

C 

C  N  on  s  + (gamma ) 

C 

DO  440,  NzMEIGHT,  NMAX,  i 

CLNM  (M,  N  )  =  PH  1 1 X  ( XH  ( 1 ,  M ) ,  XM  ( 2 ,  M ) ,  XH  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )«DELC 
440  CONTINUE 

450  CONTINUE 
C 

C  The  next  section  is  for  M  on  sd(gamma) 

C 

C  Eq  h 
C 

DO  500,  M=  MS EVN+ i ,  HEIGHT,  1 
C 

C  N  on  A (ALPHA) 

C 

DO  460,  N=  1,  MONE,  i 

IF  (DABS (XM( 1,  M)-XM( i,  N) )  . LT .  DELC/3 . 0D0 )  THEN 
CLNM (M,  N ) =CZERO 
GOTO  460 
END  IF 

CLNM  (M,  N  )  =  -  PH  1 1 X  ( XM  (  i ,  M ) ,  XM  (  2 ,  M ) ,  XM  (  i  ,  N  ) ,  XM  (  2 ,  N )  )*DELC 

460  CONTINUE 

C 

C  N  on  s- (beta  ) 

C 

DO  462,  N  =  MTW0+  1 ,  MTHR,  1 

CLNM(M,  N  )  =  C  3  *  PH  1 0  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )  *  DELC 
462  CONTINUE 

C 

C  N  on  s+ (beta  ) 

C 

DO  464,  N  =  MFOUR+  1,  MFIVE,  1 

CLNM(M,  N)  =  -C3«PHI0  (XM(  1,  M),  XM(  2,  M ) ,  XM  ( 1 ,  N  ),  XM(  2,  N )  )»DELC 
464  CONTINUE 
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c 

C  N  on  A(gamma)  and  s-(gamma) 

C 

DO  470,  N=MFIVE+ 1,  MSEVN,  i 

CLNM(M,  N )  =  PHI 1 Y  (XM(  i,  M ) ,  XM (  2,  H ) ,  XH ( 1 ,  N ),  XM(  2,  N )  )»DELC 
470  CONTINUE 

C 

C  N  on  sd(gamma) 

C 

DO  480,  N=MSEVN+ 1,  MEIGHT,  1 
CLNM(H,  N ) =CZERO 

IF  (N  .EQ.  M)  CLNM(M,  N)=  (-0.5D0,  0.D0) 

480  CONTINUE 

C 

C  N  on  s+ ( gamma ) 

C 

DO  490,  N:MEIGHT+i,  NMAX,  1 

CLNM(M,  N )  :  PH  I  i  Y  ( XM  ( 1 ,  H  ) ,  XM  ( 2,  H ) ,  XM  ( 1 ,  N  ) ,  XM  (  2,  N  )  )»DELC 
490  CONTINUE 

500  CONTINUE 
C 

C  The  last  major  section!!  Is  for  M  on  s+ (gamma) 

C 

C  Eq  i 
C 

DO  550,  M=MEIGHT+ 1,  NMAX,  i 
C 

C  N  on  A (ALPHA ) 

C 

DO  510,  N=l,  MONE,  1 

CLNM(M,  N  )  =  -  PH  1 1 Y  ( XM  ( 1 ,  M ) ,  XM  (  2,  M ) ,  XM  ( 1 ,  N  ) ,  XM(  2,  N  )  )  *  DELC 
510  CONTINUE 

C 

C  N  on  A (BETA) 

C 

DO  512,  N=MONE+  1,  MTWO,  1 

CLNM  (M,  N  )  =  -  C  3  *  PH  1 0  ( XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2,  N  )  )  »  DELC 
512  CONTINUE 

C 

C  N  on  sd (beta ) 

C 

DO  514,  N  =  MTHR+ 1 ,  MFOUR,  1 

CLNM(M,  N)=C3»PHI0(  XM  ( 1 ,  M ) ,  XM  (  2 ,  M ) ,  XM  ( 1 ,  N  ) ,  XM  (  2 ,  N  )  )*DELC 
514  CONTINUE 

C 

C  N  on  A(gamma)  and  s-(gamma) 

C 

DO  520,  N=MF I VE+ 1 ,  MSEVN,  1 

CLNM(M,  N  )  =  -  PHI  1 X  ( XM ( 1 ,  M  ) ,  XM  (  2,  M ) ,  XM (  1 ,  N  ) ,  XM  ( 2,  N  )  )«DELC 
520  CONTINUE 
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c 

C  N  on  sd(gamma) 

C 

DO  530,  N  =  MSEVN+i,  MEIGHT,  1 

CLNM  (M,  N )  =  -  PH  I  i  X  ( XM  ( 1 ,  M ) ,  XM  ( 2,  M ) ,  XM  (  i ,  N  ) ,  XH  ( 2,  N )  )»DELC 
530  CONTINUE 

C 

C  Non  s  + ( gamma ) 

C 

DO  540,  N  =  ME I GHT+ 1 ,  NMAX,  1 
CLNM (H,  N ) =CZERO 

IF  (N  .EQ.  M)  CLNM (M,  N  )  =  (  -  0 . 5D0,  0  .  DO  ) 

540  CONTINUE 

550  CONTINUE 
C 
C 

C  The  CLNM( )  matrix  is  now  filled. 

C 

C 

WRITE  («,«)  'LNM  MATRIX  FILLED. .. GOING  TO  INVERT' 

C 

C  The  major  matrices  are  now  all  filled.  The  next  step  is  to 
C  invert  the  CLNM( )  matrix. 

C 

C  Using  AFIT’s  IMSL  library,  one  statement  will  invert  CLNM(  ) 

C 

CALL  DL2NCG  (NMAX,  CLNM,  500,  CLNM,  500,  CVKS,  IWKS  ) 

WRITE  («,*)  »  MATRIX  INVERTED...' 

C 

C  We  next  fill  the  "g()"  matrix  of  known  values 
C 

DO  601,  L4=  1,  IMAX,  1 

c 

c  AL2  is  the  input  angle  in  degrees 
c 

AL2  =  DBLE( FLOAT (L4) ) *  RSTEP-RSTEP+  RFRST 
WRITE { * , 888 )  AL2 
888  FORMAT ( 5X,  F6 . 2 ) 

c 

c  Convert  the  input  angle  to  radians 
c 

THTI =  AL2*  PI/ 1 80 . DO 
c 

c  Fill  the  G ( )  matrix 
c 

DO  560,  1  =  1,  MONE,  1 

G ( I ) =CZERO 
560  CONTINUE 

DO  570,  I  =MONE+ 1 ,  MTWO,  1 
G ( I  )  = CKPO ( I , THT I ) / 2 . DO 
570  CONTINUE 
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DO  580,  I  =MTWO+  i  ,  NMAX,  1 
G  ( I  )  =  CZERO 
580  CONTINUE 

Multiply  the  inverted  matrix  by  the  input  matrix,  G() 
using  the  IMSL  library  function. 

CALL  DMCRCR  ( NMAX,  NMAX,  CLNM,  500,  NMAX,  i,  G,  500,  NMAX,  i,  COUT,  500) 
c 

c  If  it’s  the  first  input  angle,  then  write  out  the  currents 
c 

IF  (L4  .EQ.  1)  THEN 
DO  700,  IM1  :  1 ,  MONE,  1 

WRITE (* ,  860  )  REAL ( COUT ( IMi ) ),  DIMAG (COUT ( IM1 ) ) 
WRITE  (5,  860  )  REAL  ( COUT  ( IMi  )  ),  DIMAG  (COUT  ( IMI  )  ) 

860  FORMAT  (IX,  2E19.12) 

700  CONTINUE 

END  IF 

WRITE  ( 5,  *  )  »  ' 

WRITE  («,*)’  ’ 

c 

c  Calculate  the  RCS  using  Eq.  121  in  the  Thesis 
c  for  each  output  angle  from  -90  to  90 
c 

DO  600,  L  =  -  90,  90,  1 

THTR =DBLE (FLOAT (L) )*PI/180. DO 
CSUM=  0 . DO 
DO  610,  M=  1,  MONE,  1 

C SUM r  CSUM+ COUT (M ) *  CDEXP ( ~  J*  WK*  XM ( 1 ,  M ) *  DS IN ( THTR ) ) 
610  CONTINUE 

IF  ( L  . NE .  0 )  THEN 

CSUM=CSUM*  DS IN (WK*  DELC*  0 . 5D0*  DS IN (THTR ) )/ (WK* DS IN ( THTR ) ) 
END  IF 

IF  (L  .EQ.  0)  THEN 

CSUM=CSUM*  DELC/2 . DO 
ENDIF 

RCS (L4, L ) r  4 . DO*  C DABS ( CSUM ) *  CDABS (CSUM) 

600  CONTINUE 

601  CONTINUE 
c 

c  Output  the  RCS  data 
c 

DO  603,  Lr -90,  90,  1 

WRITE  (3,  8  68  )  L,  RCS(1,  L),  RCS  (2,  L),  RCS  (3,  L),  RCS  (4,  L) 
868  FORMAT  (14,  4E19.12) 

603  CONTINUE 
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c  Close  the  output  file  and  go  back  for  the  next  K  and  d 
c 

CLOSE (UNIT=  3 ) 

650  CONTINUE 

CLOSE (UNIT=2 ) 

CLOSE (UNIT  =  5 ) 

CLOSE (UNIT=4 ) 

END 


CKPO  is  the  function  that  gives  the  physical  optics  current 
across  the  aperture.  This  is  actually  the  H  field  incident 
to  the  aperture. 

COMPLEX*  1  6  FUNCTION  CKPO(M,THTI) 

IMPLICIT  DOUBLE  PRECISION  A,  B,  D,  E,  F,  H,  0,  Q,  R,  T,  U,  V,  W,  X,  Y,  Z 

IMPLICIT  COMPLEX*  16  C,  J,  G,  P 

COMMON  J,  XM(  1:  2,  1:  500  ),  WK,  DELC 

CKPO  =  2 . D0*CDEXP ( - J*WK* ( XM ( 1 ,  M ) *  DS I N ( THT I ) ) ) 

RETURN 

END 


PHIO  is  the  function  that  gives  the  Hankie  function  of  the 
first  kind  of  order  0  for  K  times  r.  It  is  multiplied  by 
j/4  to  give  the  Green’s  function  in  two  dimensions. 

COMPLEX* 16  FUNCTION  PHIO (XI,  Y1 ,  X2,  Y2 ) 

IMPLICIT  DOUBLE  PRECISION  A,  B,  D,  E,  F,  H,  0,  Q,  R,  T,  U,  V,  W,  X,  Y,  Z 

IMPLICIT  COMPLEX*  16  C,  J,  G,  P 

COMMON  J,  XM(  1:  2,  1:  500  ),  WK,  DELC 

WKR=WK*  DSQRT ( (X1-X2 )* (XI -X2 )+ (Y1-Y2 )* (Y1-Y2 ) ) 

PHIO  = (DBS JO (WKR )  + J*  DBSYO (WKR ) )* J*0 . 2  5D0 

RETURN 

END 


c  PHI1X  gives  the  x  component  of  the  gradient  of  PHIO 
c 

COMPLEX* 16  FUNCTION  PHI  1 X  ( X 1 ,  Y 1 ,  X2,  Y2  ) 

IMPLICIT  DOUBLE  PRECISION  A,  B,  D,  E,  F,  H,  0,  Q,  R,  T,  U,  V,  W,  X,  Y,  Z 

IMPLICIT  COMPLEX*  1  6  C,  J,  G,  P 

COMMON  J,  XM(  1:  2,  1:  500  ),  WK,  DELC 

R= DSQRT ( (X1-X2 )* (X1-X2 )+ (Y1-Y2 )* (Y1-Y2 ) ) 

WKR  =WK*  R 

PHI 1X  =  J*  0 . 25D0*WK* (X1-X2 )*  (DBSJ1 (WKR )  + J* DBSY1 (WKR) )/R 

RETURN 

END 
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C 

c  PHI1Y  gives  the  y  component  of  the  gradient  of  PHIO 
c 

COMPLEX* i 6  FUNCTION  PHI i Y ( X i ,  Y1 ,  X2 ,  Y2  ) 

IMPLICIT  DOUBLE  PRECISION  A,  B,  D,  E,  F,  H,  0,  Q,  R,  T,  U,  V,  W,  X,  Y,  Z 

IMPLICIT  COMPLEX*  16  C,  J,  G,  P 

COMMON  J,  XM(i:  2,  1:  500  ),  WK,  DELC 

R=DSQRT( (X1-X2 )* (X1-X2 )+ (Y1-Y2  )* (Y1-Y2  )  ) 

WKR  =WK*  R 

PH I  i  Y  =  J  *  0  . 2  5D0*WK*  (Y1-Y2  )*  ( DBS J 1 (WER )  + J* DBSY1 (VKR ) )/R 

RETURN 

END 
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Appendix  Bj  Humeri c a  1  I y  Generated  Data  For  the  Teat  Case 

The  plots  contained  in  this  appendix  were  obtained  from 
the  computer  program  in  Appendix  A.  The  program  is  discussed 
in  Chapter  7  and  some  of  the  results  are  analy2ed  in  the 
results  section,  Chapter  8. 
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Figure  17.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalised  Wave  Number,  R,  Equal  to  0.1,  a  Cavity  Depth, 
0.25,  and  an  Incident  Plane  Wave  With  0s  =  0°. 
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Figure  18.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  0.1,  a  Cavity  Depth,  d,  of 
0.25,  and  Incident  Plane  Waves  with  0j  =  22.5°,  45.0°  and 

67. 5°. C. 
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k=0.1  d— 1.0  Incident  Angle=0° 


Figure  19-  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  0.1,  a  Cavity  Depth,  d,  of 
1.0,  and  an  Incident  Plane  Wave  With  =  0°. 
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Figure  20.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normal  ized  Wave  Number,  R,  Equal  to  10.0,  a  Cavity  Depth,  d,  of 
1.0,  and  an  Incident  Plane  Wave  With  0j  =  0°. 
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Figure  21.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  0.1,  a  Cavity  Depth,  d,  of 
1.0,  and  Incident  Plane  Waves  with  6j  =  22.5°,  45.0°  and  67.5°. 
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Figure  22.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  O.i,  a  Cavity  Depth,  d,  of 
4.0,  and  an  Incident  Plane  Wave  With  0j  =  0°. 
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Figure  23.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  R,  Equal  to  0.1,  a  Cavity  Depth,  d,  of 
4.0.  and  Incident  Plane  Waves  with  0i  =  22.5°.  45.0°  and  67.5° 
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Figure  24.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  O.i,  a  Cavity  Depth,  d,  of 
8.0,  and  an  Incident  Plane  Wave  With  9^  =  0°. 
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Figure  25.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  Jc,  Equal  to  0.1,  a  Cavity  Depth,  d,  of 
8.0,  and  Incident  Plane  Waves  with  0*  :  22.5°,  45.0°  and  67.5°. 
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Figure  26.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  1 . 0,  a  Cavity  Depth,  d,  of 
0.25,  and  an  Incident  Plane  Wave  With  =  0°. 
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Figure  27.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  R,  Equal  to  i.O,  a  Cavity  Depth,  d,  of 
0.25,  and  Incident  Plane  Waves  with  0j  =  22.5°,  45.0°  and 

67.5°. 
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Figure  28.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  i.O,  a  Cavity  Depth,  d,  of 
1.0,  and  an  Incident  Plane  Wave  With  0j  =  0°. 
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Figure  29.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  1.0,  a  Cavity  Depth,  d,  of 
8.0,  and  an  Incident  Plane  Wave  With  0j  =  0°. 


121 


k=  1 .0  d=4.0  Incident  Angle=0° 


Figure  30.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  i.O,  a  Cavity  Depth,  d,  of 
4.0,  and  an  Incident  Plane  Wave  With  9j  =  0°. 


122 


k  <r 


I 

I 

I 

I 

I 

I 

I 

i 


k=1.0  d=4.0 


Bistatic  Scattering  Angle  (degrees) 


Figure  31.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normal  ized  Wave  Number,  K,  Equal  to  1.0,  a  Cavity  Depth,  d,  of 
4.0,  and  Incident  Plane  Waves  with  0^  =  22.5°,  45.0°  and  67.5°. 
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Figure  32.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  10.0,  a  Cavity  Depth,  d, 
0.25,  and  an  Incident  Plane  Wave  With  0j  =  0°. 
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Figure  33.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Norma  1  i  zed  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth,  d,  of 
0.25,  and  an  Incident  Plane  Wave  With  =  22.5°. 
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k=10.0  d=0.25  Incident  Angle=45° 


Figure  34.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth,  d,  of 
0.25,  and  an  Incident  Plane  Wave  With  0^  =  45°. 
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Figure  35.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  10.0,  a  Cavity  Depth,  d, 
0.25,  and  an  Incident  Plane  Wave  With  =  22.5°,  45°,  and 

67.5°. 
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Figure  36.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  k,  Equal  to  10.0,  a  Cavity  Depth,  d, 
4.0,  and  an  Incident  Plane  Wave  With  0j  =  22.5°,  45°,  and 

67.5°. 
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Figure  37.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth, 
4.0,  and  an  Incident  Plane  Wave  With  eA  =  0°. 
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Figure  30.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth,  d, 
4.0,  and  an  Incident  Plane  Wave  With  0j  =  22.5°. 
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k=10.0  d=4.0  Incident  Angle=45° 


Figure  39.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth, 
4.0,  and  an  Incident  Plane  Wave  With  0i  =  45.0°. 
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Figure  40.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth,  d,  of 
1.0,  and  an  Incident  Plane  Wave  With  8j  :  22.5°,  45°,  and 

67.5°. 
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igure  42.  Bistatic  RCS  of  a  Rectangular 
Normalized  Wave  Number,  K,  Equal  to  10.0, 
10,  and  an  Incident  Plane  Wave  With  6^  = 
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Figure  43.  Bistatic  RCS  of  a  Rectangular  Cavity  With  a 
Normalized  Wave  Number,  K,  Equal  to  10.0,  a  Cavity  Depth,  d, 
1.0,  and  an  Incident  Plane  Wave  With  0*  =  45.0  . 
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Bistatic  RCS  at  0°/45 
is  demonstrated  again 
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Figure  48.  Bistatic  RCS  at  0°/45°  with  Varying  Depth  and  K=6.0 
Reciprocity  is  demonstrated  once  again. 
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Appendix  C_:  A  Model  for  Obtaining  Measurements 

A  scatterer  has  been  designed  for  a  physical  problem 
related  to  the  example  problem  attacked  in  this  thesis. 

Perhaps  a  follow-on  thesis  could  concentrate  on  the 
measurements . 

Measuring  the  fields  scattered  from  the  channel  is  not  a 
trivial  matter.  The  plane  on  which  the  incident  fields  impinge 
has  been  assumed,  thus  far,  to  be  infinite.  Such  an  infinite 
plane  that  is  also  perfectly  conducting  cannot  be  constructed. 
However,  some  approximations  can  be  made  to  at  least  simplify 
those  problems.  Further,  proper  calibrations  while  making  the 
measurements  will  also  help  create  accurate  measurements. 

If  the  plane  containing  the  channel  to  be  measured  cannot 
be  infinite,  then  it  should  be  very  large;  but  how  large?  If 
the  edges  of  the  plane  are  in  the  far  field,  then  it  should  be 
sufficiently  large  enough  for  the  incident  wave  to  "see"  a 
cavity  in  an  "infinite"  plane.  The  far  field  is  loosely 
defined  as  a  distance,  r,  such  that  (22:21) 

2nr  >>  X 

For  most  practical  applications,  the  far  field  begins  at  about 
1 0 X  from  the  scatterer  (9).  In  this  problem,  the  scatterer 
is  the  cavity.  If  the  edges  of  the  plane  were  at  least  10X 


from  the  cavity  center,  they  may  be  considered  to  be  infinite. 


The  assumption  of  infinite  conductivity  needs  to  be 
considered  when  making  measurements .  Though  perfect  conductors 
do  exist  (i.e.  superconductors),  it  is  not  yet  practical  to  use 
them  in  most  applications.  A  good  approximation  to  a  perfect 
conductor  for  problems  of  this  type  would  be  any  "good" 
conductor.  Aluminum  is  a  good  conductor  (li).  For  measurement 
purposes  it  can  be  considered  to  be  a  PEC  (9).  A  plane  and 
cavity  made  of  aluminum  can  therefore  be  approximated  as 
perfect  conductors. 

When  making  the  measurements,  the  system  needs  to  be 
calibrated  to  nullify  anomalies  that  may  exist  within  the 
system  (such  as  transmitter/receiver  peculiarities,  mounting 
scatter,  etc.).  If  the  system  were  calibrated  using  the  above 
plane  without  the  cavity,  then  most  of  the  edge  effects  of  the 
plane  itself  can  be  subtracted  from  the  measurements  along  with 
the  system  anomalies.  This  will  mostly  leave  only  measurements 
of  the  fields  due  to  the  cavity  alone.  This  is  the  same 
quantity,  the  integral  in  Eq  (76),  that  has  been  calculated  and 
plotted  in  the  earlier  sections  of  this  thesis. 

Using  the  described  target  and  the  described  technique,  it 
should  be  possible  to  duplicate  the  data  plotted  in  the  last 
section.  The  target  has  been  designed  and  built  and  is  the 
property  of  the  Air  Force  Institute  of  Technology;  to  date,  no 
measurements  have  been  made. 

In  addition  to  a  very  large  aluminum  plane,  an  adjustable 
channel  width  and  depth  is  also  desired.  Using  all  of  the 
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